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Tensegrity  is  an  abbreviation  of  tension  and  integrity.  Tensegrity  structures  are 
spatial  structures  fonned  by  a  combination  of  rigid  elements  in  compression  (struts)  and 
connecting  elements  that  are  in  tension  (ties).  In  three-dimensional  tensegrity  structures 
no  pair  of  struts  touches,  and  the  end  of  each  strut  is  connected  to  non-coplanar  ties, 
which  are  in  tension.  In  two-dimensional  tensegrity  structures,  struts  still  do  not  touch.  A 
tensegrity  structure  stands  by  itself  in  its  equilibrium  position  and  maintains  its  fonn 
solely  because  of  the  arrangement  of  its  struts  and  ties.  The  potential  energy  of  the 
system  stored  in  the  springs  is  at  a  minimum  in  the  equilibrium  position  when  no  external 
force  or  torque  is  applied.  A  closed-form  solution  of  a  two-spring,  three-spring,  and  four- 
spring  planar  tensegrity  mechanism  was  developed  to  determine  all  possible  equilibrium 
configurations  when  no  external  force  or  moment  is  applied.  Here  “closed  form”  means 
that  all  solution  equilibrium  poses  will  be  determined,  although  for  each  case  a  high- 
degree  polynomial  will  have  to  be  solved  numerically. 
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CHAPTER  1 
LITERATURE  REVIEW 


We  examined  the  literature  related  to  tensegrity  systems,  self-deployable  tensegrity 
systems,  and  parallel  mechanisms  with  compliant  elements. 

Knight  et  al.  [1]  reported  on  the  line  geometries  of  a  family  of  tensegrity  structures 
called  skew  prisms  (anti-prisms,  tensegrity  prisms)  with  pairs  of  triangles,  squares, 
pentagonals,  hexagonals,  and  so  on,  located  at  tops  and  bases.  They  used  the  quality 
index  with  S-P-S  connectors  (S  is  a  ball-and-socket  joint,  and  P  is  a  sliding  joint).  The 
quality  index  measures  the  geometric  stability  of  in-parallel  devices  (which  in  three 
dimensions  consist  of  a  pair  of  rigid  platforms  connected  by  legs  that  are  kinematical  S- 
P-S  connectors).  Geometric  stability  depends  on  the  geometry  of  the  legs. 

Figures  1-1  A  and  1-1  C  show  a  sequence  of  parallel  prisms  with  parallel  ties  of 
fixed  length. 


(A) 


(C) 


Figure  1-1.  Parallel  prisms. 
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A  relative  anti-clockwise  rotation  can  be  achieved  by  rotating  the  top  of  each  prism 
through  its  vertical  axes  as  shown  in  Figures  1-2  B  and  1-2  D  to  yield  a  corresponding 
right-handed  tensegrity  prism  that  is  completed  by  inserting  struts  on  the  diagonals  of  the 
skew  quadrilaterals.  Left-handed  tensegrity  prisms  are  mirror  images  of  right-handed 
prisms  and  are  obtained  by  rotating  the  tops  in  a  clockwise  direction  and  interchanging 


the  ties  and  struts.  The  angle  a  is  unique  for  each  tensegrity  prism  (a  =  90 


180 


) 


n 

where  n  is  the  total  number  of  sides  in  the  upper  or  lower  polygons.  The  value  of  a 
(together  with  the  size  of  the  tops  and  bases  and  their  distance  apart)  enables  one  to 
compute  the  length  of  the  struts  and  ties  that  defines  a  unique  configuration  for  the 
tensegrity  prism. 


(B) 


Figure  1-2.  Tensegrity  prisms. 

The  quality  index  for  each  skew  prism  is  zero,  which  means  that  it  has 
instantaneous  mobility.  Further  examination  revealed  that  all  the  sets  of  connector-lines 
for  each  prism  belong  to  a  linear  complex  of  lines  within  a  5 -system  of  screws.  Each  of 
the  sets  of  lines  is  reciprocal  to  a  single  screw.  Adding  additional  ties  along  the  diagonals 
of  the  skew  prism  faces  can  easily  form  reinforced  skew  prisms. 
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Unloaded  tensegrity  prisms  that  are  stable  in  the  sense  that  they  are  configurations 
of  minimum  potential  energy,  but  they  have  instantaneous  mobility  simply  because  the 
connector-lines  belong  to  a  linear  complex.  Reinforced  tensegrity  prisms  where 
additional  ties  are  inserted  are  completely  stable  and  shown  in  Figure  1-3. 


Tensegrity  Prisms  Reinforced  Tensegrity  Prisms 


Figure  1-3.  Plan  view  of  tensegrity  and  corresponding  reinforced  prisms. 

Pellegrino  [2]  shows  how  deployable  structures  are  able  to  change  their  shape  from 
a  packed  configuration  to  operational  form.  Usually,  energy  is  stored  in  the  structure 
when  packed  and  released  during  unpacking  when  the  operational  configuration  is 
required.  Some  simple  examples  are  a  spring  loaded  umbrella,  a  retractable  roof  of  a  car, 
and  a  space  radio  telescopic  antenna.  Deployable  structures  are  used  for  ease  of 
transportation  and  storage.  The  essential  requirement  is  that  the  transformation  process 
should  be  autonomous  and  reliable,  and  without  causing  any  damage  to  nearby  structures. 

Examples  of  deployable  structures  that  undergo  large  geometric  transfonnations  are 
coiled  rods,  flexible  shells,  lattice  column,  and  membranes.  A  deployable  structure  using 
lattice  columns  with  several  longitudinal  elements,  called  longerons,  braced  at  regular 
intervals  by  short  members  perpendicular  to  the  longerons  and  by  diagonal  members,  has 
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low  dead  weight  but  can  support  high  load.  This  structure  has  a  very  small  wind 
resistance  in  atmosphere  and  minimal  particle  meteoroid  damage  in  space. 

Tarnai  [3]  presents  a  brief  analogy  between  compatibility  and  equilibrium  of  a 
finite  linkage  mechanism.  Some  linkages  are  known  whose  degree  of  freedom  in  some 
specific  position  is  greater  than  the  expected  degree  of  freedom  as  calculated  from  the 
equilibrium  equation.  An  example  is  a  four  bar  linkage  with  equal  length  opposite  bars. 
The  linkage  has  two  shapes:  one  associated  with  a  parallelogram  and  the  other  with  an 
anti  parallelogram.  Plotting  these  angles  as  a  function  of  each  other  provides  two  curves 
with  a  common  point.  This  point  is  when  all  four  bars  he  on  a  straight  line.  At  this 
instant  the  mechanism  may  move  from  the  parallelogram  configuration  to  the  anti¬ 
parallelogram  configuration. 

Duffy  et  al.  [4-6]  analyze  a  three  dimensional  tensegrity  structure  that  is  made  up  of 
elastic  and  rigid  elements  (Figure  1-4).  In  this  assembly,  the  elastic  elements  are  under 
tension  and  the  rigid  elements  under  compression.  The  papers  present  the  static  position 
analysis  problem  and  detennine  the  position  assumed  by  the  structure  when  external 
loads  are  applied,  and  when  the  system  is  presented  by  changing  the  free  lengths  of  the 
compliant  elements. 
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Figure  1-4.  Component  of  a  tensegrity  structures. 

The  mathematical  formulation  to  find  the  equilibrium  positions  of  the  structure  is 
based  on  the  virtual  work  principle.  The  obtained  equations  are  solved  using  numerical 
methods.  Some  assumptions  are  made  to  simplify  the  derivation  of  the  mathematical 
model.  These  assumptions  are  the  absence  of  internal  dissipative  forces  and  the  manner 
by  which  the  external  forces  are  applied.  The  numerical  method  determines  the 
coordinates  of  the  strut  end  points  in  the  equilibrium  position.  A  force  balance  is  then 
conducted  to  validate  the  results. 

Stern  [7]  presents  the  position  analysis  of  a  symmetric  n-strut  tensegrity  system.  A 
three  dimensional  n-struts  system  with  two  platforms,  one  at  the  top  and  the  other  at  the 
bottom,  is  considered  in  the  analysis.  A  static  analysis  of  the  internal  forces  is  conducted 
on  the  top  and  bottom  platforms.  The  relationship  between  the  geometry  of  the  structure 
and  the  internal  forces  are  investigated.  Tensegrity  structures  with  different  number  of 
struts  are  analyzed.  The  results  of  each  analysis  are  compared  with  the  results  of  other 
analyses  to  obtain  common  patterns  in  all  systems  by  relating  the  results  to  the  number  of 
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struts  in  the  systems.  The  patterns  are  formulated  into  equations  based  on  the  number  of 
struts  in  the  system. 

The  formulation  is  done  for  3,  4,  5,  and  6  struts  tensegrity  systems  with  solid  top 
and  bottom  platforms.  For  example,  a  plan  view  of  4  strut  tensegrity  system  is  shown  in 
Figure  1-5. 

ID!  GG 


Figure  1-5.  4-strut  tensegrity,  3D,  top  view. 

Marshall  [8]  introduces  a  parallel  platform  device  that  incorporates  tensegrity 
principles.  The  device,  shown  in  Figure  1-5,  replaces  the  struts,  of  a  tensegrity  structure 
with  prismatic  actuators,  and  each  elastic  member  with  a  cable-spring  combination  in 
series.  The  length  of  the  three  prismatic  actuators  and  the  length  of  the  cables  that  are  in 
series  with  the  spring  are  adjustable  and  thus  the  device  has  six  degrees  of  freedoms.  This 
study  shows  that  in  order  to  achieve  an  arbitrary  position  and  orientation  of  the  top 
platform,  an  external  wrench  must  be  applied  to  maintaining  equilibrium.  This  study  also 
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shows  that  the  device’s  compliance  characteristics  can  be  varied  while  maintaining  its 
position  and  orientation. 

A  reverse  analysis  of  the  device  is  presented  in  which  the  desired  pose  and  total 
potential  energy  are  given  and  the  lengths  of  the  prismatic  actuators  and  the  cables  are 
determined.  The  effect  of  a  seventh  leg,  another  prismatic  actuator,  is  also  analyzed  and 
found  to  satisfactorily  implement  the  needed  external  wrench  (Figure  1-6). 


Figure  1-6.  Tensegrity  platform. 

Gantes  and  Konitopoulou  [9]  discusses  bi-stable  deployable  structures.  A  bi-stable 
structure  is  self  standing  and  stress  free  when  fully  closed  or  fully  deployed.  It  exhibits 
incompatibilities  between  the  member  lengths  at  intermediate  geometric  configurations 
during  the  deployment  process,  which  leads  to  the  occurrence  of  second-order  strains  and 
stresses  resulting  in  a  snap-through  phenomenon  that  "locks"  the  structures  in  their 
deployed  configuration.  Until  now  the  geometric  shapes  that  were  possible  in  the 
deployed  configuration  were  only  flat  or  curved  with  constant  curvature.  This  limitation 
is  addressed  in  the  paper  by  proposing  a  geometric  design  methodology  for  deployable 
arches  of  arbitrary  curvature,  accounting  also  for  the  discrete  joint  size,  and  applying  it 
successfully  for  the  geometric  design  of  a  semi-elliptical  arch.  The  arch  is  then  modeled 
with  finite  elements,  and  a  geometrically  nonlinear  analysis  is  performed  in  order  to 
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verify  the  deploy  ability  feature.  Further  verification  is  provided  by  the  construction  of  a 
small-scale  physical  model.  A  preliminary  structural  design  indicates  the  overall 
feasibility  of  the  arch  for  short  to  medium  spans  and  light  loads. 

Deployable  structures  are  prefabricated  space  frames  consisting  of  straight  bars 
linked  together  in  the  factory  as  a  compact  bundle,  which  can  then  be  unfolded  into  large- 
span,  load  bearing  structural  shapes  by  simple  articulation.  Because  of  this  feature  they 
offer  significant  advantages  in  comparison  to  conventional,  non-deployable  structures  for 
a  wide  spectrum  of  applications  ranging  from  temporary  structures  to  the  aerospace 
industry,  being  mainly  characterized  by  their  feature  of  transforming  and  adapting  to 
changing  needs. 

From  a  structural  point  of  view,  deployable  structures  have  to  be  designed  for  two 
completely  different  loading  conditions,  under  service  loads  in  the  deployed 
configuration,  and  during  deployment.  The  structural  design  process  is  very  complicated 
and  requires  successive  iterations  to  achieve  some  balance  between  desired  flexibility 
during  deployment  and  desired  stiffness  in  the  deployed  configuration  From  a  geometric 
point  of  view,  the  whole  idea  of  this  type  of  deployable  structure  is  based  on  the  so-called 
scissor-like  elements,  pairs  of  bars  connected  to  each  other  at  an  intennediate  point 
through  a  pivotal  connection  which  allows  them  to  rotate  freely  about  an  axis 
perpendicular  to  their  common  plane  but  restrains  all  other  degrees  of  freedom,  while,  at 
the  same  time,  their  end  points  are  hinged  to  the  end  points  of  other  scissor-like  elements. 

Yin  et  al.  [10]  present  a  special  planar  three-spring  mechanism  that  is  designed  to 
control  contact  forces.  An  energy  function  is  defined  to  describe  the  behavior  of  this  kind 
of  mechanism  and  it  can  be  used  to  perfonn  the  catastrophe  analysis  of  this  mechanism. 
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The  analysis  result  can  be  used  as  a  design  and  control  tool.  By  comparing  the  three- 
spring  system  and  a  two-spring  system,  it  was  found  that  the  three-spring  mechanism  has 
better  stability  than  the  two-spring  system.  A  three-spring  mechanism  which  can  be  used 
to  control  a  general  contact  force  in  a  plane  is  also  analyzed. 

Intuitively,  Yin  et  al.  [10]  showed  a  catastrophe  occurs  whenever  a  smooth  change 
of  parameters  gives  rise  to  a  discontinuous  change  in  behavior.  A  well  known  example 
that  can  easily  be  made  to  demonstrate  a  catastrophe  is  Zeeman’s  catastrophe  machine 
(see  Figurel-7).  Zeeman’s  machine  can  be  constructed  by  attaching  two  linear  springs  to 
a  single  point  C  on  a  disk  that  can  rotate  about  0.  One  of  the  springs  is  attached  to  a  fixed 
pivot  at  point  A  and  the  other  spring  is  attached  to  point  B  which  can  be  moved  in  the  x, 
y  plane.  The  position  of  point  B  is  called  the  controlling  parameter  as  it  dictates  the 
position  of  the  disk  which  is  defined  by  the  angle  between  OA  and  OM. 


Figure  1-7.  Zeeman’s  catastrophe  machine. 

Parallel  compliant  mechanisms  are  classically  known  for  their  role  in  mounting  and 
suspension  systems.  Recently,  however,  they  have  been  instrumental  in  the  field  of 
robotic  force  control.  A  new  theory  for  the  simultaneous  control  of  displacement  and 
force  for  a  partially  constrained  end-effecter  has  been  proposed  in  Yin  et  al.  [10], 
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Figure  1-8  shows  a  planar  compliant  mechanism  which  is  actually  a  planar  two- 
spring  system.  It  consists  of  a  pair  of  linear  springs  connected  at  one  end  to  a  movable 
base  and  at  the  other  end  to  a  common  pivot  which  is  the  axis  of  the  wheel  contacting  to  a 
surface.  The  base  is  connected  to  a  planar  two  freedom  P-P  (P  denotes  prismatic  pair  ) 
manipulator  .  The  contact  force  can  be  controlled  by  displacing  the  two  prismatic  joints 
of  the  manipulator.  The  required  displacements  can  be  calculated  from  the  stiffness 
mapping.  This  kind  of  control  was  called  kinestatic  control  by  Griffis  and  Duffy  [1 1].  In 
order  to  design  the  planar  two-spring  system  it  is  necessary  to  compute  a  spring  stiffness 
which  will  generate  a  range  of  displacements  of  the  movable  base  which  can  be  produced 
by  the  prismatic  joints  over  a  required  range  of  change  in  contact  force.  Clearly  if  the 
system  is  over-designed  and  the  spring  stiffness  is  very  high  it  is  always  possible  to 
generate  any  necessary  changes  in  contact  force.  However  such  a  system  will  be  too 
sensitive  to  errors  because  very  small  displacements  of  the  platform  will  generate  large 
changes  in  contact  force.  On  the  other  hand  if  the  springs  are  too  soft  there  can  be 
stability  problems. 


Figure  1-8.  Planar  two-spring  system. 

Yin  et  al.  [10]  are  also  considered  another  compliant  mechanism.  That  is  shown  in 
Figure  1-9.  It  consists  of  three  linear  springs  joined  to  the  triangular  frame  fixed  points 
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and  connected  at  the  axis  of  the  wheel.  The  triangular  frame  is  connected  to  the  planar 
two  freedom  P-P  manipulator.  This  mechanism  is  a  special  three-spring  system.  The 
catastrophe  analysis  of  this  system  was  presented.  Comparing  the  results  demonstrates 
that  three-spring  system  has  better  stability  characteristics  than  the  two-spring  system. 


Figure  1-9.  Special  three-spring  system. 

A  more  general  planar  three-spring  compliant  mechanism  is  shown  in  Figure  1-10. 
This  mechanism  is  connected  to  the  planar  three  freedom  R-R-R  (  R  denotes  revolute 
pair)  manipulator.  It  can  be  used  to  control  both  force  and  moment.  The  catastrophe 
analysis  of  this  mechanism  was  perfonned  in  Yin  et  al.  [10]. 


*x/7 

/  NY 


Figure  1-10.  General  planar  three-spring  compliant  mechanism. 

Crane  and  Duffy  [12]  studied  kinematics  of  robot  manipulators  and  Crane  and 
Bayat  [13]  presented  cases  of  two  dimensional  tensegrity  structures  that  proved  to  be 
mathematically  challenging  in  spite  of  their  simple  configurations.  Different  methods 
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were  tried  to  solve  these  problems.  It  was  expected  the  final  closed  form  solution  be  a 
simple  mathematical  expression.  In  each  case  the  intermediate  mathematical  expressions 
become  complicated  to  solve  and  the  relevant  matrices  become  large  to  manipulate. 

Roth  and  Whiteley  [14]  propose  a  technology  based  on  tensegrity  for  tough,  rigid, 
large-scale  domes  that  are  also  economical  to  construct.  The  development  of  a  structural 
technology  to  economically  cover  large  areas  would  be  useful  for  warehouses,  permanent 
or  temporary  protection  for  archaeological  and  other  vulnerable  sites,  large-scale 
electrical  or  electromagnetic  shielding  and  exclusion  or  containment  of  flying  animals  or 
other  objects.  Structures  based  on  such  a  technology  can  serve  as  frameworks  in  which 
environmental  control,  energy  transformation  and  food  production  facilities  could  be 
embedded.  The  space  application  is  also  possible  by  using  self-deployed  structures. 
Summary  Advantages  are  improved  rigidity,  ethereal,  resilient,  equal-length  struts, 
simple  Joints. 

The  resulting  list  represents  an  initial  attempt  to  identify  applications  for  which  the 
technology  would  be  suitable.  As  the  technology  develops  and  is  tested  against  them, 
some  or  all  of  the  applications  may  be  winnowed  out  and  other  suitable  applications  not 
in  the  list  may  become  apparent.  The  current  list  is  as  follows: 

Superstructures  for  embedded  substructures  allow  the  substructures  to  escape 
terrestrial  confines  where  this  is  useful  (e.g.  in  congested  or  dangerous  areas,  urban  areas, 
flood  plains  or  irregular,  delicate  or  rugged  terrains). 

1 .  Economic  large-scale  protection  of  storage,  archaeological,  agricultural, 
construction  or  other  sites. 

2.  Refugee  or  hiking  shelters. 

3.  Frames  over  cities  for  environmental  control,  energy  transfonnation  and  food 
production. 
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4.  Large-scale  electrical  or  electromagnetic  shielding. 

5.  Exclusion  or  containment  of  flying  animals  or  other  objects. 

6.  Spherical  superstructures  for  space  stations. 

7.  Earthquake-resistant  applications. 

These  structures  are  extremely  resilient  and  testing  would  very  likely  show  they 
could  withstand  large  structural  shocks  like  earthquakes.  Thus,  they  would  likely  be 
desirable  in  areas  where  earthquakes  are  a  problem.  A  tensegrities,  pneumatic,  structure 
performed  very  well  during  recent  earthquakes  in  Japan. 

•  Low-environmental-impact  shells  for  musical  performances. 

•  Indoor/outdoor  pavilions  for  trade  shows  etc. 

•  Supports  to  hold  sunscreen  protection  for  vulnerable  amphibians. 

•  Watersheds  to  keep  rain  water  from  percolating  through  contaminated  soils  into 
groundwater,  perhaps  temporarily  during  in-situ  remediation. 

•  Frames  for  hanging  plants  or  other  objects  to  dry. 

•  Pergola,  trellis,  or  topiary  framework. 

•  Micro-meteorite  protection,  sun-shielding  for  Martian  colonies. 


Figure  1-11.  A  representation  of  a  dome  which  utilitizes  tensegrity  solutions'  technology 
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Skelton  et  al.  [15]  present  a  solution  so  the  tensegrity  structures  are  reduced  to 
linear  algebra  problems,  after  first  characterizing  the  problem  in  a  vector  space  where 
direction  cosines  are  not  needed.  That  is,  we  describe  the  components  of  all  members  by 
vectors  as  opposed  to  the  usual  practice  of  characterizing  it  as  static  problem  in  terms  of 
the  magnitude  of  tension  and  compression.  While  this  approach  transfer  the  problem  into 
vector  space  to  describe  the  problem,  the  advantage  is  that  the  vector  space  makes  the 
mathematical  structure  of  the  problem  amenable  to  linear  algebra  treatment  for  both  two 
and  three  dimensional  tensegrity  structures. 

This  paper  characterizes  the  static  equilibrium  some  tensegrity  structures. 
Furthermore,  it  uses  vectors  to  describe  each  element  (springs  and  studs),  eliminating  the 
need  to  use  direction  cosines  and  the  subsequent  transcendental  functions  that  followthe 
use  of  trigonometry  equations. 

In  this  paper,  the  authors  choose  to  represent  a  tensegrity  structure  as  an  oriented 
graph  in  real  three  dimensional  space  R3  defined  in  tenns  of  nodes  and  directed  branches 
which  are  all  represented  as  vectors  in  R3.  A  loop  is  any  closed  path  in  the  graph. 

The  advantage  of  this  approach  is  that  the  both  the  magnitude  and  the  direction 
cosines  of  the  forces  are  contained  in  vectors  which  can  be  solved  using  linear  algebra. 
Thus  linear  algebra  plays  a  larger  role  in  this  approach  compared  to  the  usual  approach  in 
mechanics  and  finite  element  methods  using  direction  cosines. 

In  this  oriented  graph,  the  nodes  consist  of  the  ends  of  the  bars.  Hence  if  there  are  n 
bars,  then  there  are  2n  nodes.  There  are  two  types  of  directed  branches;  the  string  branch 
(in  vectors)  and  the  bar  branches  (in  vectors). 
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Geometric  connectivity  of  the  structure  maintained  such  that  each  directed  branch 
can  undergo  a  displacement  in  reaching  its  equilibrium  state.  String  vectors  can  change 
both  their  length  and  orientation  while  bar  vectors  can  only  change  their  orientation. 

Node  vectors  can  change  both  their  length  and  orientation  but  subject  to  a  Law  of 
Geometric  Connectivity  is  stated  as  follows:  The  vector  sum  of  all  branch  vectors  in  any 
loop  is  zero.  These  loop  equations  are  in  the  form  of  a  set  of  linear  algebraic  equations  in 
the  branch  vectors. 

In  this  study  of  tensegrity  structures,  force  equilibrium  are  such  that  spring  can  only 
take  tension  and  bars  sustain  compressive  forces.  We  therefore  choose  to  distinguish 
between  the  string  (or  tensile)  forces  and  the  bar  (or  compressive)  forces  which  are 
defined  in  tenns  of  the  string  and  bar  vectors  respectively. 

This  paper  reduces  the  study  of  the  tensegrity  equilibrium  to  a  series  of  linear 
algebra  problems  using  directed  graph  theory.  Of  course  the  existence  conditions  for  the 
linear  algebra  problems  are  nonlinear  in  the  design  variables.  The  presented  procedure 
and  formulation  give  some  insight  to  solution  of  tensegrity  structures  and  identifies  the 
free  parameters  that  may  be  used  to  achieve  desired  structural  shapes. 

Tibert  [16]  presents  adjustable  tensegrity  structure.  A  tensegrity  structure  is  a 
lightweight  consisting  of  compression  members  surrounded  by  a  network  of  tension 
members.  They  can  be  easily  dismantled  providing  possibilities  for  reusable  and  modular 
structures.  Tensegrities  adapt  their  shape  by  self  adjusting  their  tension  and  compression 
in  their  members.  They  can  adapt  to  changing  environment  when  they  are  equipped  with 


sensors  and  actuators. 
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A  full-scale  prototype  of  an  adjustable  tensegrity  is  built  and  tested  at  Swiss 
Federal  Institute  of  Technology.  This  paper  begins  with  a  description  of  important 
aspects  of  the  design,  assembly,  and  static  testing.  Tests  show  that  the  structure  behaves 
linearly  when  subjected  to  vertical  loads  applied  to  a  single  joint.  Nonlinearities  are 
detected  for  small  displacements,  for  loads  applied  to  several  joints  and  for  adjusting 
combinations  of  telescoping  compression  members.  To  predict  behavior,  dynamic 
relaxation — a  nonlinear  method — has  been  found  to  be  reliable.  Appropriate  strut 
adjustments  found  by  a  stochastic  search  algorithm  are  identified  for  the  control  goal  of 
constant  roof  slope  and  for  the  load  conditions  studied.  When  adjusting  struts,  an 
excessive  number  of  adjustable  members  does  not  necessarily  lead  to  improved 
performance. 

Tensegrity  is  an  abbreviation  of  tensile-integrity  as  termed  by  one  R.  Buckminster 
Fuller.  He  described  tensigrity  as  “small  islands  of  compression  in  a  sea  of  tension.” 
This  is  a  description  of  a  network  of  light  tension  members  that  provide  rigidity  to  a 
limited  number  of  discontinuous  compression  members.  Some  researchers  have  found 
that  discontinuous  compressions  are  not  a  necessity  for  creating  a  tensegrity  and  more 
efficient  structures  can  be  built  if  compression  elements  are  allowed  to  join.  Defining  a 
tensegrity  structure  as  any  self-stressing  structure  composed  of  struts  and  cables  would 
include  structures  such  as  a  bicycle  wheel.  In  order  to  add  precision  tenninology,  a 
definition  has  been  proposed  by  Motro  and  Raducanu.  They  describe  a  tensegrity  as  a 
stable  system  that  contains  a  discontinuous  set  of  components  in  compression  inside  a 
network  of  components  in  tension. 
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Several  independent  efforts  have  contributed  to  the  invention  and  development  of 
tensegrity  structures.  Le  Ricolais  designed  many  unusual  lightweight  cable-strut 
structures.  The  sculptor  Kenneth  Snelson  started  in  1948  with  sculptures  employing 
members  in  continuous  tension  and  discontinuous  compression.  Applications  of  the 
tensegrity  principle  in  nature  have  recently  been  found  through  studies  that  use  a 
tensegrity  model  to  describe  observations  of  how  cells  respond  to  stress. 

Transforming  tensegrity  from  sculptures  into  practical  structural  has  been  a 
challenge  that  began  with  several  studies  of  their  geometric,  nonlinear  behavior  and  states 
of  self  stress.  Reasons  for  the  lack  of  test  data  are  related  to  the  fabrication  assembly 
process  as  well  as  to  the  lack  of  rigidity  unless  pre-stressed  members  are  used. 

Joint  design  is  the  biggest  challenge  in  constructing  a  full-scale  tensegrity. 

Although  it  is  the  problem  for  all  space  structures,  tensegrities  present  particular 
challenges.  Joints  need  to  be  pin  jointed,  modular,  and  light  in  order  to  take  advantage  of 
tensegrities  ease  of  dismantling  and  potential  reuse. 

In  addition  to  the  complexities  of  constructing  the  structure  itself,  simulation  of  the 
behavior  is  not  straightforward.  Tensegrity  structures  exhibit  geometric  nonlinear 
behavior.  In  addition,  as  it  is  common  with  other  pin-jointed  structures,  nodal  friction 
should  be  included  in  the  simulation  of  a  truss  structure  for  precision  control.  Finally,  for 
full-scale  construction,  assembly  sequences,  fabrication  tolerances,  and  construction  lack- 
of-fit  need  to  be  considered. 

Fest  et  al.  [17]  present  a  study  of  adjustable  tensegrity  structures.  A  tensegrity 
structure  is  consisting  of  compression  members  surrounded  by  a  network  of  tension 
members.  A  tensegrity  structure  is  a  lightweight  structure.  Tensegrity  structures  can  be 
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easily  dismantled  providing  possibilities  for  reusable  and  modular  structures.  Tensegrities 
adapt  their  shape  by  self  adjusting  their  tension  and  compression  in  their  members.  They 
can  adapt  to  changing  environment  when  they  are  equipped  with  sensors  and  actuators. 

A  full-scale  prototype  of  an  adjustable  tensegrity  is  built  and  tested  at  Swiss 
Federal  Institute  of  Technology.  This  paper  presents  important  aspects  of  the  design, 
assembly,  and  static  testing.  Tests  show  that  the  structure  behaves  linearly  when  it  is 
subjected  to  vertical  loads  applied  to  a  single  joint.  Nonlinearities  are  observed  when 
loads  applied  to  several  joints  and  for  adjusting  combinations  of  telescoping  compression 
members.  Appropriate  strut  adjustments  found  by  a  search  algorithm  for  the  load 
conditions  studied.  When  adjusting  struts,  an  excessive  number  of  adjustable  members  do 
not  necessarily  lead  to  equilibrium. 

Active  tensegrity  structures  have  the  potential  to  widen  the  scope  for  innovative, 
lightweight,  and  reusable  structural  systems.  Lessons  learn  through  the  construction  and 
testing  of  an  adjustable  full-scale  prototype  can  be  used  in  construction  and  design  of 
practical  tensegrity  structures.  This  work  expected  to  contribute  to  the  development  of 
tensigrity  and  improve  in  their  design  and  perfonnance  during  service. 


CHAPTER  2 

TWO  SPRING  PLANAR  TENSEGRITY  SYSTEM 

2.1  Introduction 

This  Chapter  presents  two  approaches  to  solve  the  forward  position  analysis  of  a 
two-strut  tensegrity  system  with  two  compliant  ties  and  two  non-compliant  ties  (Figure  2- 
1). 


Figure  2-1.  Two  dimensional  tensegrity  structure. 

Here  it  is  assumed  that  the  lengths  of  the  struts  and  noncompliant  ties  are  known 
together  with  the  spring  constants  and  free  lengths  of  the  two  compliant  ties.  The 
objective  is  to  determine  all  possible  equilibrium  poses  for  the  device  when  no  external 
loads  are  applied.  Gravity  loads  are  neglected. 

The  objective  of  this  effort  is  to  determine,  in  closed-form,  all  possible  equilibrium 
deployed  positions  in  stable  condition  (minimum  potential  energy)  of  a  planar  tensegrity 
system  wherein  two  of  the  ties  are  compliant.  Figure  2-2  shows  the  system  which  is 
comprised  of  two  struts  (compression  members  a n  and  a.34),  two  non-compliant  ties 
(tension  members  &4\  and  a23),  and  two  elastic  tensile  members  (springs),  one  connected 
between  points  1  and  3  and  one  between  points  2  and  4.  It  should  be  noted  in  Figure  2-2 
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that  strut  a 34  passes  through  a  slit  cut  in  strut  a!2  and  as  such  the  two  struts  do  not 
intersect  or  collide. 


Figure  2-2.  Two-spring  planar  tensegrity  system. 

2.2  Problem  Statement 

The  problem  statement  for  the  two-spring  plane  tensegrity  system  can  be  explicitly 
written  as 
given: 

ai2,  a.34  lengths  of  struts 

a23,  a4i  lengths  of  non-compliant  ties 

ki,  L01  spring  constant  and  free  length  of  compliant  tie  between  points  4  and  2 
k2,  L02  spring  constant  and  free  length  of  compliant  tie  between  points  3  and  1 
find: 

Li  length  of  spring  1  at  equilibrium  position 
L2  length  of  spring  2  at  equilibrium  position 

It  should  be  noted  that  the  problem  statement  could  be  formulated  in  a  variety  of 
ways,  that  is,  a  different  variable  (such  as  the  relative  angle  between  strut  a34  and  tie  a4! ) 
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could  have  been  selected  as  the  generalized  parameter  for  this  problem.  Two  solution 
approaches  are  presented. 

2.3  Approach  1:  Determine  (Li,  L2)  to  Minimize  Potential  Energy 
2.3.1  Development  of  Geometric  Equations 

Figure  2-3  shows  the  three  angles  04,  04,  and  0M4  which  must  satisfy  the  relation 


04  +  041  —  n  +  04M  . 


(2-1) 


Figure  2-4  shows  the  triangle  formed  by  side  a34,  a23,  and  Li.  A  cosine  law  for  this 
triangle  can  be  written  as 


2 


(2-2) 


3 


a23 


Figure  2-3.  Identification  of  angles  04,  04,  and  0M4- 
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Figure  2-4.  Triangle  4-1-3. 
Solving  for  cos04'  yields 


COS04'  = 


2  _ i  2  _  2 

a23  ~^1  ~a34 

2  L  i  a34 


(2-3) 


Figure  2-5  shows  the  triangle  formed  with  sides  a4i,  an,  and  Li.  A  cosine  law  for 
this  triangle  can  be  written  as 


+^—  +  Li  a41  cos04"  = 


2  2 
Solving  for  cos04"  yields 


(2-4) 


2  _  T  2  _  2 

OH  ai2  ~~  —  a41 

4  = 


(2-5) 


2Lj  a41 

Figure  2-6  shows  the  triangle  formed  by  sides  a4i,  a34  and  L2.  A  cosine  law  for  this 
triangle  can  be  written  as 


22  j  : 

^34  ^41  r\  ^  2 

~  + - +  a34  a41  cos  04  =  — — - 


2  2 
Solving  for  cos04  yields 


T  2  2  2 

n  L 2  ^34  ^41 

cos  =  — - — - — 


(2-6) 


2  ^34  d/\\ 


(2-7) 
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Figure  2-5.  Triangles  4-1-2. 


Figure  2-6.  Triangles  4-1-3. 
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Equating  the  cosine  of  the  left  and  right  sides  of  2-1  yields 

cos  (04  +  04')  =  cos  (n  +  04")  .  (2-8) 

Expanding  this  equation  yields 

cos04  cos04'  -  sin04  sin04’  =  -  cos04"  .  (2-9) 

Rearranging  2-9  yields 

cos04  cos04'  +  cos04"  =  sin04  sin04’  (2-10) 

Squaring  both  sides  of  2-10  gives 

(cos04)2  (COS04’)2  +  2  cos04  cos04'  cos04"  +  (COS04")2  =  (sin04)2  (sin04')2  .  (2-11) 

Substituting  for  (sin04)  and  (sin04')  in  terms  of  cos04  and  cos04’  gives 
(cos04)2  (COS04’)2  +  2  COS04  COS04'  COS04"  +  (COS04")2  = 

(l-cos204)(l-cos204’)  .  (2-12) 

Equations  2-3,  2-5,  and  2-7  are  substituted  into  2-12  to  yield  a  single  equation  in 
the  parameters  Li  and  L2  which  can  be  written  as 

AL24  +  BL22  +  C  =  0  (2-13) 

where 

A  =  Li2,  B  =  Li4  +  B2Li2  +  B0,  C  =  C2Li2  +  C0  (2-14) 

and 

B2  =  -  (a232  +  a34“  +  a4i2  +  ai22), 

Bo  =  (ai2  -  a4i)  (ai2  +  a4i)  (a23  -  a34)  (a23  +  a34)  , 

C2  =  (a34  -  a4i)  (a34  +  a4J)  (a23  -  ai2)  (a23  +  ai2)  , 

Co  =  (a4ia23  +  a34ai2)  (a4Ja23-  a34ai2)  (a4i  +  a23"  -  ai2~  -  a34~) 


(2-15) 
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Equation  2-13  expresses  length  L2  as  a  function  of  the  length  Li.  The  function  is 
quadratic  with  respect  to  (L2)2  and  (Li)2  and  thus  there  will  be  two  possible  values  for 
(L2)2  for  each  value  of  (Li)2. 

As  a  verification  of  Equation  2-13  a  numerical  example  is  presented  to  show  how 
four  values  of  L2  can  occur  for  a  given  value  of  Li.  For  this  example,  assume  a4J  =  1.0  m, 
ai2  =  2.0  m,  a23  =  2.0  m,  and  a34  =  2.5  m.  Further,  assume  the  value  of  Li  is  given  as  Li  = 
2.25  m.  Evaluating  the  coefficients  in  Equations  2-15  and  2-14  and  substituting  into 
Equation  2-13  gives 

(5.0625)  L24  +  (-58.3242)  L22  +  (1 10.2500)  =  0  .  (2-16) 

Solving  this  equation  for  L2  yields  four  answers; 

L2a  = -3.0228,  L2b  =  3.0228,  L2c  = -1.5438,  L2d=  1.5438  .  (2-17) 

Figure  2-7  shows  the  mechanism  in  four  configurations  where  for  each  of  these 
configurations  L2  can  have  a  positive  or  negative  value,  thereby  giving  eight  possible 
states  of  the  mechanism.  Flowever,  due  to  the  symmetry  associated  with  the  reflected 
solutions  about  a4i,  only  four  values  of  L2  (or  two  values  of  L2  )  exist.  This  example  was 
presented  to  show  that  the  degree  of  the  polynomial  which  relates  Li  and  L2  is  indeed 


fourth  order. 
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Figure  2-7.  Possible  configurations  for  numerical  example. 

The  potential  energy  of  the  system  can  be  evaluated  as 

U  =  tk,(L,-L01)2+ik2(L2-Lm)2 

2  2  .  (2-18) 

At  equilibrium,  the  potential  energy  will  be  a  minimum.  This  condition  can  be 

determined  as  the  configuration  of  the  mechanism  whereby  the  derivative  of  the  potential 

energy  taken  with  respect  to  the  length  Li  equals  zero,  i.e. 


dU 

dL, 


—  kj  (Lj  L01)  +k2  (L,  L02) 


dL, 

dLj 


(2-19) 


The  derivative  dlAdLi  can  be  determined  via  implicit  differentiation  from 
Equation  2-13  as 


dL2 

dL, 


-L,[L22(L22 


r,  T  2  2  2  2 

2Lj  a  23  a41  a34  a12 


2)  +  (a122  - a2 


:  Xa4l  a3 


)] 


L2[Lj2  (L,2 


TT  2  2  2  2 

'  2L2  a  23  a41  a34  a12 


2)  +  (a122  -a4 


,2)(a  232  -a342))] 


(2-20) 
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Substituting  Equation  2-20  into  Equation  2-19  and  regrouping  gives 

D  L25  +  E  L24  +  F  L23  +  G  L2~  +  H  L2  +  J  =  0  (2-21) 

where 

D  =  Dj  Li,  E  =  E,  Li,  F  =  F3  Li3  +  F2  Li2  +  Fi  U,  G  =  G3  Li3  +  Gi  U, 

H  =  H5  L!5  +  H4  L!4  +  H3  Lj3  +  H2  L!2  +  Hj  U  +  H0,  J  =  Ji  U  (2-22) 

and  where, 

Di  =  k2,  Ei  =  -k2  L02, 

F3  =  2  (k2  -ki),  F2  =  2  ki  L01,  Fi  =  -  k2  (ai22  +  a232  +  a342  +  a4i2), 

G3  =  -2  k2  L02,  Gi  =  k2  L02  (an  +  a23“  +  a34  +  a4f), 

H5  =  -  ki,  H4  =  ki  L01,  H3  =  ki  (ai2“  +  a232  +  a34“  +  a4f ), 

H2  =  -  ki  L01  (ai2“  +  a332  +  a342  +  a4f ), 

Hi  =  -  kl  (a342  -  a232)  (a4i2  -  ai22)  +  k2  (a342  -  a4i2)  (a232  -  ai22), 

Ho  =  ki  L01  (a34  -  a23")  (a4f  -  ai2~), 

Ji  =  k2  L02  (a34“  -  a4E)  (ai2“  -  a332)  (2-23) 

Equations  2-13  and  2-21  represent  two  equations  in  the  two  variables  Li  and  L2. 

The  simultaneous  solution  of  these  equations  is  a  necessary  condition  that  the  system  is  in 
equilibrium. 

2.3.2  Sylvester’s  Solution  Method 

A  system  of  m  equations  of  order  n  in  the  variable  Y  may  be  written  as  (Gantes  and 
Konitopoulou  [9]) 

Ai,n  Yn  +  Ai,n.i  Yn_1  +  Ai,n_2  Yn'2  +  . . .  +  Ai,o  =  0 
A2,n  Yn  +  A2,n_,  Yn_1  +  A2,„-2  Yn'2  +  . . .  +  A2,o  =  0  , 

Am,n  Yn  +  Am>n.i  Yn-‘  +  Am,n.2Yn'2  +  . . .  +  Am,o  =  0  .  (2-24) 
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The  terms  Un  =  Yn  ,  Un-i  =  Yn_l  , .  .  . ,  Ui  =  Y,  Uo  =  1  may  be  substituted  into  Equation 
2-24  to  yield 


Al,n  Un  +  Al>n-1  Un-1  +  Aijn.2  Un-2  +  •  • 

.  +  A IX)  Uo  =  0 

A2,n  Un  +  A2jn-1  Un-1  +  A2,n-2  Un-2  +  •  • 

•  +  A2,o  Uo  =  0  , 

Amjn  Un  +  Am,n-1  Un-1  +  Amin-2  Un_2  +  . 

•  •  +  Am  0  Uo  =  0  . 

(2-25) 

Equation  2-25  represents  a  system  of  m  homogeneous  equations  with  n+1  unknowns  (the 
tenn  Uo  will  be  treated  as  an  unknown  for  the  time  being)  and  can  be  written  in  matrix 
form  as 
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If  the  number  of  equations  is  less  than  the  number  of  unknowns,  additional 
equations  can  be  created  by  multiplying  selected  equations  from  set  of  Equation  2-24  by 
powers  of  Y  to  obtain  a  new  system.  For  example,  multiplying  the  first  equation  from  the 
set  of  Equation  2-24  by  Y  will  introduce  one  new  equation  and  one  new  unknown,  Yn+1 
(which  will  be  represented  by  Un+i).  Now,  multiplying  the  second  equation  of  the  set  of 
Equation  2-24  by  Y  will  yield  another  equation,  but  no  new  variables  are  introduced. 
Ultimately,  the  system  will  be  formed  where  the  new  total  number  of  homogeneous 
equations  equals  the  new  total  number  of  unknowns.  In  other  words,  the  coefficient 
matrix  will  be  a  square  matrix.  Now,  one  solution  to  the  set  of  ‘homogeneous’  equations 
is  the  trivial  solution  where  all  the  U;  terms  equal  zero.  However,  this  cannot  occur  since 
Uo=l .  Thus  it  must  be  the  case  that  the  set  of  equations  are  linearly  dependent  which 
means  that  the  determinant  of  the  coefficient  matrix  must  equal  zero. 
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Now  consider  that  the  terms  Ay  of  the  coefficient  matrix  are  polynomial  functions 
of  order  q  of  another  variable  called  X,  that  is 

Ai  =  a;  Xq  +  bi  Xq-‘  +  Ci  xq'2  +  ...  +  Zi  =  0  (2-27) 

where  a; ,  b; ,  c; ,  . . .  z;  are  constants.  Evaluation  of  the  determinant  of  the 
coefficient  matrix  will  result  in  a  polynomial  in  the  variable  X  that  must  equal  zero  so 
that  the  set  of  equations  in  the  variables  U;  will  be  linearly  dependent.  This  solution 
approach  will  now  be  applied  to  Equations  2-13  and  2-21  which  represent  two  equations 
in  the  two  unknowns  Li  and  L2. 

2.3.3  Solution  of  Geometry  and  Energy  Equations 

Equations  2-13  and  2-21  can  be  solved  by  using  Sylvester’s  variable  elimination 

2  3  4 

procedure  by  multiplying  Equation  2-13  by  L2,  L2  ,  L2  ,  and  L2  and  Equation  2-21  by  L2, 

2  3* 

L2",  L2  to  yield  a  total  of  nine  equations  that  can  be  written  in  matrix  form  as 
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This  set  of  equations  can  only  be  solved  if  the  determinant  of  the  9X9  coefficient 
matrix  is  equal  to  zero.  Expansion  of  this  detenninant  yields  a  30th  degree  polynomial  in 
the  variable  Li.  When  the  determinant  was  expanded  symbolically,  it  was  seen  that  the 
two  lowest  order  coefficients  were  identically  zero.  Thus  the  polynomial  can  be  divided 
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throughout  by  Li  to  yield  a  28th  degree  polynomial.  The  coefficients  of  the  28th  degree 
polynomial  were  obtained  symbolically  in  terms  of  the  given  quantities.  They  are  not 
presented  here  due  to  their  length  and  complexity. 

Values  for  L2  that  correspond  to  each  value  of  Li  can  be  detennined  by  first  solving 
Equation  2-13  for  four  possible  values  of  L2.  Only  one  of  these  four  values  also  satisfies 
Equation  2-21. 

2.3.4  Numerical  Examples 

Example  1 :  The  following  parameters  were  selected  for  this  numerical  example: 

Strut  lengths  ai2  =  3  in.  a34  =  3.5  in. 

Non-compliant  tie  lengths  a.41  =  4  in.  a23  =  2  in. 

Spring  1  free  length  &  spring  constant  L0i  =  0.5  in.  ki  =  4  lbf/in. 

Spring  2  free  length  &  spring  constant  L02  =  1  in.  k2  =  2.5  lbf/in. 

Eight  real  and  twenty  complex  roots  were  obtained  for  the  28th  degree  polynomial 
in  Li.  Real  values  for  Li  and  the  corresponding  values  of  L2  are  shown  in  Table  2-1. 

Table  2- 1 .  Eight  real  solutions. 

Case  1  2  3  4  5  6  7  8 

LI,  in.  -5.485  -5.322  -1.741  -1.576  1.628  1.863  5.129  5.476 

L2,  in.  2.333  -2.901  -1.495  1.870  1.709  -1.354  -3.288  2.394 

The  values  of  Li  and  L2  listed  in  Table  2-1  satisfy  the  geometric  constraints  defined 
by  Equation  2-13  and  the  energy  condition  defined  by  Equation  2-21.  Each  of  these  eight 
cases  was  analyzed  to  determine  whether  it  represented  a  minimum  or  maximum 
potential  energy  condition  and  cases  3,  4,  5,  and  6  were  found  to  be  minimum  states.  A 
free  body  analysis  of  struts  ai2  and  a34  was  performed  to  show  that  these  bodies  were 
indeed  in  equilibrium  for  each  of  these  four  cases.  Figure  2-8  shows  the  four  static 
equilibrium  configurations  of  the  system.  Note  that  there  are  several  cases  where  the 
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Case3 


2 


spring  in  compression  with 
a  negative  spring  length 


spring  in  tension 


Figure  2-8.  Four  equilibrium  configurations. 

spring  lengths  are  negative,  such  as  for  example  the  spring  between  points  2  and  4  in  case 
4.  In  this  case  the  spring  is  acting  to  pull  point  2  towards  point  4  which  is  counter 
intuitive  for  the  nonnal  spring  that  is  in  compression. 

There  are  twenty  complex  roots  for  Li  in  this  example.  Unique  corresponding 
values  for  L2  were  obtained  for  each  case  (Table  2-2).  In  fact  there  are  corresponding 
values  for  L2  such  that  the  complex  pair  of  Li  and  L2  satisfies  both  Equations  2-13  and 
2-2 1 .  This  result  means  that  no  extraneous  tenns  were  introduced  in  the  elimination 


procedure.  These  complex  roots  were  not  analyzed  further. 


Table  2-2.  Complex 

solutions. 

Case 

Li,  inches 

L2,  inches 

1 

-7.0285  -0.0648i 

-0.3225  +  2.5486i 

2 

-7.0285  +  0.0648i 

-0.3224 -2.5486i 

3 

-2.0919 -2.4639i 

-6.5995  +  0.3617i 

4 

-2.0919 +  2.4639i 

-6.5995  -  0.36 1 7i 

5 

-1.7904  -2.0044i 

6.5232  +  0.0738i 

6 

-1.7904 +  2.0044i 

6.5232 -0.0738i 

7 

-0.7211  -0.1827i 

0.2411  -  1.2909i 

8 

-0.7211  +  0.1827i 

0.2411  +  1.2909i 

9 

0.005823  -  1.7651i 

-0.05803  -0.01682i 

10 

0.005823  +  1.7651i 

-0.05803  +  0.01682i 

11 

0.3062 -0.410H 

-0.1220 -0.9493i 

12 

0.3062 +  0.410H 

-0.1220 +  0.9493i 

13 

0.8140 -0.5553i 

0.5055  +0.954H 
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Table  2-2.  Continued. 


Case 

Li,  inches 

L2,  inches 

14 

0.8140  +  0.5553i 

0.5055  +0.954H 

15 

2.3264  -  1.9808i 

6.1894  +  0.2416i 

16 

2.3264+  1.9808i 

6.1894 -0.2416i 

17 

2.9910 -2.5413i 

-6.2069- 0.93  lOi 

18 

2.9910 +  2.5413i 

-6.2069  + 0.93  lOi 

19 

7.0351  -0.1067i 

-0.4518 -2.6442i 

20 

7.0351  +  0.1067i 

-0.4518 +  2.6442i 

Example  2:  For  this  case,  a  C  language  program  was  written  that  would  randomly 
select  values  for  the  parameters  an,  a23,  a34,  a4i,  ki,  L0i,  k2,  and  L02  in  attempt  to  find  a 
case  where  the  number  of  real  roots  for  Li  was  twenty  eight.  Such  a  case  was  not  found, 
but  a  set  of  inputs  resulting  in  twenty  four  real  values  for  Li  is  presented  here. 

The  input  parameters  for  this  case  were  selected  as: 

Strut  lengths  a^  =  3.309848  in.  a34  =  3.002692  in. 

Non-compliant  tie  lengths  a4i  =  7.335484  in.  a33  =  7.978210  in. 

Spring  1  free  length  &  spring  constant  L0i  =  0.953703  in.  ki  =  0.193487  lbf/in. 

Spring  2  free  length  &  spring  constant  L02  =  0.6073 18  in.  k2  =  9. 108249  lbf/in. 

Twenty  four  real  and  four  complex  roots  were  obtained  for  Li.  The  real  values  for 
Li  and  the  corresponding  values  of  L2  are  shown  in  Table  2-3.  The  four  complex  values 


for  Li  were  1.2818  ±  18.2760i  and  -0.8019  ±  18.2940L 
Table  2-3.  Twenty-four  real  solutions. _ 


Case 

Li,  inches 

L2,  inches 

1 

-19.0226 

-2.5540 

2 

-16.9651 

2.8637 

3 

-11.7045 

0.6215 

4 

-11.0691 

4.3323 

5 

-11.0649 

-4.3325 

6 

-10.5742 

-4.6685 
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Table  2-3.  Continued. 


Case 

Li,  inches 

L2,  inches 

7 

-10.5716 

4.6687 

8 

-5.1348 

10.3382 

9 

-5.1343 

-10.3382 

10 

-3.7712 

-11.2881 

11 

-3.7704 

11.2881 

12 

-0.02315 

11.7686 

13 

-0.02315 

-11.7686 

14 

3.7733 

11.2881 

15 

3.7737 

-11.2881 

16 

5.1330 

-10.3382 

17 

5.1334 

10.3382 

18 

10.5735 

4.6686 

19 

10.5755 

-4.6685 

20 

11.0631 

-4.3326 

21 

11.0664 

4.3324 

22 

11.7046 

0.6193 

23 

17.4121 

2.7902 

24 

19.5242 

-2.4884 

At  first  glance,  it  appears  that  a  set  of  input  parameters  have  been  found  whereby 
most  of  the  roots  of  the  resulting  polynomial  equation  in  Li  and  the  corresponding  values 
of  L2  are  real.  However,  from  Figure  2-2,  it  is  apparent  that  the  system  will  be  realizable 
only  if 

|a4i  -  ai2|  ^  |Li|  <  a4i  +  a[2  , 

|a34  -  a23|  <  |Li|  <  a34  +  a23 , 

|a23  -  a42|  <  |L2|  <  a23  +  ai2  , 


|a4i  -  a34|  <  |L2|  <  a4i  +  a34  . 


(2-29) 
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For  this  numerical  case,  the  system  will  be  realizable  only  if 
4.9755  in.  <  |Li|  <  10.6453  in.  , 

4.6684  in.  <  |L2|  <  10.3382  in.  (2-30) 


Table  2-4  lists  the  cases  which  satisfy  these  conditions. 


Table  2-4.  Eight  feasible  real  solutions. 


Case 

Li,  inches 

L2,  inches 

6 

-10.5742 

-4.6685 

7 

-10.5716 

4.6687 

8 

-5.1348 

10.3382 

9 

-5.1343 

-10.3382 

16 

5.1330 

-10.3382 

17 

5.1334 

10.3382 

18 

10.5735 

4.6686 

19 

10.5755 

-4.6685 

The  spring  L2  is  at  an  extreme  limit  value  in  every  one  of  the  cases  listed  in  Table 
2-4.  When  |L2|  =  10.3382  in.,  strut  a34,  tie  a4i,  and  the  spring  L2  are  collinear.  When  |L2| 
=  4.668,  strut  a[2,  tie  a23,  and  spring  L2  are  collinear.  Both  configurations  are  shown  in 
Figure  2-9. 


2 


Figure  2-9.  Realizable  configurations  for  case  2. 
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2.3.5  Discussion  of  Results 

It  is  apparent  that  obtaining  values  of  Li  and  L2  that  simultaneously  satisfy 
Equations  2-13  and  2-21  is  a  necessary  condition  for  equilibrium.  Satisfying  these 
equations  is  not  sufficient,  however,  to  guarantee  that  the  system  is  physically  realizable. 

The  real  values  for  Li  and  L2  that  were  calculated  in  the  second  example,  but  which 
violate  the  conditions  of  Equation  2-29  are  an  interesting  case.  Here  some  or  all  the 
angles  84,  84'  and  84"  will  be  complex,  yet  the  condition  of  Equation  2-1  can  still  be 
satisfied. 

Presented  above  was  the  technique  to  obtain  all  possible  equilibrium  positions  of  a 
planar  tensegrity  system  that  incorporates  two  compliant  members.  The  approach  of 
satisfying  geometric  constraints  while  simultaneously  finding  positions  where  the 
derivative  of  the  total  potential  energy  with  respect  to  the  generalized  coordinate  equaled 
zero  resulted  in  a  28th  degree  polynomial  in  a  single  variable.  Although  the  resulting 
polynomial  was  of  higher  degree  than  anticipated,  an  analysis  of  the  real  and  complex 
solutions  indicates  that  no  extraneous  solutions  were  introduced  during  the  variable 
elimination  procedure.  Complex  solutions  satisfy  the  geometry  equation  but  they  can  not 
be  used  to  construct  the  structure,  hence  they  are  not  the  solution. 

2.4  Approach  2:  Determine  (cos  84  and  cos  0j)  to  Minimize  Potential  Energy 

The  objective  of  this  approach  is  to  investigate,  in  closed-form,  possible 
equilibrium  positions  using  the  cosine  of  two  angles,  84  and  81,  as  the  descriptive 
parameters  for  the  system. 

Figure  2-18  shows  the  tensegrity  system  which  is  comprised  of  two  struts 
(compression  members  ai2  and  a34),  two  non-compliant  ties  (tension  members  a4i  and 
a23),  and  two  elastic  tensile  members  (springs),  one  connected  between  points  1  and  3  and 
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one  between  points  2  and  4.  It  should  be  noted  in  Figure  2-10  that  strut  a 34  passes 
through  a  slit  cut  in  strut  ai2  and  as  such  the  two  struts  do  not  intersect  or  collide. 


Figure  2-10.  Tensegrity  system  with  two  struts,  two  ties  and  two  springs. 
The  problem  statement  is  written  as 
given  a4i,  ai2,  a23,  a34 


ki,  k2,  L01,  Lo2 

find  cos04  (and  corresponding  value  of  cosOi)  when  the  system  is  in 
equilibrium 

The  solution  approach  is  outlined  as  follows: 

1 .  Obtain  expression  for  cosOi  in  terms  of  cos04 

2.  Write  Li  in  terms  of  cos0i  and  L2  in  terms  of  cos04 

3.  Write  the  potential  energy  equation 

4.  Obtain  an  equation  in  the  variables  cos04  and  cos0i  which  corresponds  to  the 
condition  dU/dcos04  =  0 

5.  Utilize  Sylvester’  method  to  obtain  values  for  cos04  and  cos0i  that  satisfy  the 
equation  from  step  1  and  step  4. 

2.4.1  Obtain  Expression  for  cos0i  in  terms  of  cos04 

The  cosine  law  for  the  quadrilateral  1 -2-3-4  is  written  as 
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7  —  a23 

^41  -  ~ 


where  by  definition 


(2-31) 


■Z4i  —  a12  (X4S[  +  Y4Cj  )  +  Z4  + 


X4  —  a34  s4 

Y4  =  -  (a4i  +  a34  c4) 


Z4=^-  +  ^-  +  a34a41c4 


Substituting  Equation  2-3 1  into  2-32  and  rearranging  gives 


ai2  a  23 


(2-32) 

(2-33) 

(2-34) 

(2-35) 


(2-36) 


ai2X4Si  —  a12Y4Cj  +Z4  + 

Substituting  for  X4,  Y4,  and  Z4  and  then  squaring  this  equation  and  substituting  s4z=  1  -c4z 
2  2 

and  Si  =1-Ci  yields 


A  Ci  +  B  Ci  +  D  -  0 


(2-37) 


where 


A  —  Ai  c4  +  A2 

2 

Ai  =  -2  ai2  a34  a4i 

2  2  2  2 
A2  =  -ai2“  a34“  -  a  12“  a4i 


B  —  Bj  c4“  +  B2  c4  +  B3 
2 

Bi  =  -  2  an  a4i  a34 

2  2  2  2 
Bo  =  ai2  a34  (-a3Y  -  ai2~  -  3  a4i  +  ao3  ) 

2  2  2  2 
B3  =  an  a4i  (-a4i  -  ai2~  -  a34"  +  a23~) 
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D  -  Di  C4  +  D2  C4  +  D3 

Di  -  -  a34  a4i  -  ai2  a34 

D2  =  a34  a4i  (-a4f  -  ai22  -  a34"  +  a232) 

D3  =  (-ai22  +  a23~  —  2  ai2  a34  —  a34~  —  a4p)  (-ai22  +  a23~  +  2  ai2  a34  —  a34~  —  a4i2)  /  4 


(2.38) 


Equation  (2-37)  provides  a  relationship  between  ci  and  C4. 


2.4.2:  Write  Li  in  terms  of  cosOi  and  L2  in  terms  of  cos04 


A  cosine  law  for  triangle  4-1-2  that  expresses  Li  in  terms  of  ci  may  be  written  as 


^  +  ^A  +  a12a4lCl=^ 
2  2  12  41  1  2 


(2-39) 


A  cosine  law  for  the  triangle  3-4-1  that  expresses  L2  in  terms  of  C4  may  be  written  as 


^34  ,  &41  ,  ^7 

2  2  34a4lC4  =^“ 


(2-40) 


2.4.3:  Write  the  potential  energy  equation 

The  energy  stored  in  the  springs  is  given  by 


u  =  y2  ki  (Li-Loi)2  +  y2  k2  (L2-L02)2  ■ 


(2-41) 


2.4.4:  Express  the  equation  dU/dcos04  =  0  in  terms  of  the  variables  c4  and  c4 

The  derivative  of  the  potential  energy  U  with  respect  to  c4  may  be  written  as 


dLj  dcj 
dc,  dc4 


d  -k2(L2-L02)‘ 


(2-42) 


^  =  k,(L,  -  Lm)  Al  A  +  kj(L2  -  l02)  Al 
dc4  dct  dc4  dc4 


(2-43) 


From  Equations  2-39  and  2-40 
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dL[  _  a12  a4i 
dct  Lj 


(2-44) 


dL2  _  &34  ^4i 

dc4  L2 
From  Equation  2-37 

dCj  _  Aj  Cj2  +  2Bj  Cj  c4  +  B2  Cj  +  2Dj  c4  +  D2 
dc4  2  Aj  Cj  c4  +  2  At  Cj  +  Bj  c4~  +  B2  c4  +  B3 

Substituting  Equations  2-44,  2-45,  and  2-46  into  2-43  gives 


(2-45) 


(2-46) 


dU-WT  I  \ 

—  kj(Lj  L01) 
dc, 


f  \ 

ai2a41 

Aj  Cj  +  2Bj  Cj  C^  +  B2  Cj  4"  2Dj  C^j  +  D2 

+k2(L2  L02) 

f  \ 

a34a41 

V  J 

^2  Aj  Cj  C4  +  2A2  Cj  +  Bj  C^j  +  B2  C4  +  B3  J 

l  ^2  J 

(2-47) 


Equating  Equation  2-47  to  zero,  dividing  throughout  by  a4i  and  rearranging  gives 


kj(Lj  L01)L2  a12F  +  k2(L2  L02)Lja34G  0 


(2-48) 


where 


F  —  Aj  Cj  +2 Bj  Cj  c4  +  B2  Cj  +  2D)  c4  +  D2 
G  =  2  Aj  Cj  c4  +  2  A,  Cj  +  Bj  c4~  +  Bt  c4  +  B3 
Rearranging  this  Equation  2-48  gives 


(2-49) 

(2-50) 


Lj  L2  [kj  a12  F  +  k2  a34  G]  —  Lj  [k2  L02  a34  G]  +  L,  [kj  L01  a)2  F] 
Squaring  both  sides  and  rearranging  gives 


(2-51) 


Lj2  L22  [kj  a12  F  +  k,  a34  G]2  -Lj2  [k2  L0,  a34  G]2  -L22  [kj  L0,  a12  F]2  = 


2Lj  L2  [k2  L02  a34  G][kj  L01  a12  F] 
Squaring  both  sides  again  gives 


(2-52) 
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L,4  L24  [k,  a12  F  +  k2  a34  G]4  +  L,4  [k2  L02  a34  G]4  +  L24  [k,  L01  a12  F]4 
-  2 Lj4  L22  [k,  a12  F]2  [k2  L02  a34  G]2  -  2  L,2  L24  [k,  a12  F]2  [k,  L01  a12  F]2 
+  2L/  L22  [k2  L02  a34  G]2  [k,  L01  a12  F]2 
=  4L12L22[k2L02a34G]2[k1L01a12F]2 


(2-53) 


Rearranging  gives, 


Li4  L24  [k,  a12  F  +  k2  a34  G]4  +  L,4  [k2  L02  a34  G]4  +  L24  [k,  L01  a12  F]4 
-  2 Lj4  L22  [k,  a12  F]2  [k2  L02  a34  G]2  -  2  L,2  L24  [k,  a12  F]2  [k,  L01  a12  F]2 
_2Lj2  L22  [k2  L02  a34  G]2  [k,  L01  a12  F]2  =  0 


(2-54) 


2  o 

Substituting  for  Li  and  L2“  using  Equations  2-39  and  2-40  and  factoring  the 
polynomial  in  terms  of  C4  and  Ci  gives 


C10  ci10  +  C9  ci9  +  C8  ci8  +  C7  ci7  +  C6  ci6  +  C5  ci5 

+  C4  ci4  +  C3  ci3  +  C2  ci2  +  Cici  +  Co  =0  (2-55) 


where 


C10  -  Cio,2  c4  +  C  10,1  C4  +  Cio,o 

C9  =  C9;3  C43  +  C9j2  C4“  +  C9,i  C4  +  C9,o 

C8  =  Cgi4  C44  +  Cs,3  C43  +  Cs,2  C4“  +  Csj  C4  +  Cs,0 

C7  =  C7,5  C45  +  C7,4  c44  +  C7,3  c43  +  C7,2  c42  +  C7,1  c4  +  C7.0 

C6  =  C6,6  C46  +  C6,5  C45  +  C6,4  C44  +  C6,3  C43  +  C6,2  C4“  +  C6,l  c4  +  C6,0 

C5  =  C5J  C47  +  €5,6  C46  +  €5,5  C45  +  C5j4  C44  +  C53  C43  +  C5,2  C4-  +  C5,1  C4  +  €5,0 

C4  =  C4,8  c48  +  C4,7  c47  +  C4,6  c46  +  C4j5  c45  +  C4j4  c44  +  C4,3  c43  +  C4,2  C4”  +  C4,i  C4  + 


C4,o 


C3  -  C3,9  C49  +  C3;8  C48  +  C3j  C47  +  C3j6  C46  +  C3  5  C45  +  C3j4  C44  +  C3j3  C43  +  C3,2  c4“  + 

c3,l  c4  +  C3,o 
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C2  —  C2,10  C410  +  C2,9  C49  +  C2,8  C48  +  C2,7  C47  +  C2,6  C46  +  C2,5  C4'  +  C2,4  C £  +  C2,3  C 4  3  + 
C2,2  C4”  +  C2,l  C4  +  C2,0 

Cl  =Ci,io  c410  +  Cl  9  C49  +  Ci,8  C48  +  Cl  ,7  C4  +  Cl, 6  C46  +  Cl, 5  C4~  +  Cl, 4  C44  +  Cl, 3  C43  + 


C  1,2  C4”  +  Ci,i  C4  +  Ci,o 

Co  =  Co, 10  C410  +  Co, 9  C49  +  Co, 8  C48  +  Co, 7  C47  +  Co, 6  C46  +  Co, 5  C45  +  Co,4  C44  +  C0,3  C43  + 
Co, 2  C42  +  Co,i  C4  +  Co  (2-56) 

All  the  terms  Cy  have  been  obtained  symbolically  and  are  defined  in  terms  of 
known  mechanism  parameters.  For  example, 


Cio,o=  8  ki4  ai26  a4i4  a342Ai4  +  4ki4  ai26  a4i6  A[4+  4ki4  ai26  a4i2  a344  Ai4 
Cio.i  =  16ki4  ai26  a4i5  a34  Ai4  +  16  ki4  ai26  a4J3  a343  Ai4 

Ci 0,2  =  16  ki4  an6  a4i4  a342  Ai4  c42  (2-57) 


The  remaining  terms  are  quite  lengthy  and  are  not  listed  here. 

2.4.5  Construct  Sylvester  Matrix  for  this  approach 

23456789 

Multiplying  Equation  2-37  by  Ci,  Ci  ,  ci  ,  ci  ,  ci  ,  Ci  ,  ci  ,  Ci  ,  ci  and  2-55  by  ci 
gives  12  equations  that  can  be  written  as  follows 


"  0 

0 

0 

0 

0 

0 

0 

0 

0 

A 

B 

D 

11 

Ci 

"0" 

0 

0 

0 

0 

0 
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A 
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D 

0 

10 

Cl 

0 

0 

0 

0 

0 

0 

0 

0 

A 
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D 

0 

0 
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C, 

0 

0 

0 

0 

0 
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C, 
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0 
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0 
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0 

0 

0 

0 

4 

Cl 

0 

0 

A 

B 

D 

0 

0 

0 

0 

0 

0 

0 

0 

3 

Cl 

0 

A 

B 

D 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

Cl 

0 

0 

Cio 

c9 

c8 

c7 

C6 

C5 

C4 

C3 

C2 

Ci 

Co 

Cl 

0 

c 

c9 

Cg 

c7 

C6 

C5 

C4 

C3 

C2 

C, 

Co 

0 

1 

0 

.  (2-58) 
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Expansion  of  the  12x12  determinant  yields  a  32nd  degree  polynomial  in  the  parameter  C4. 
2.4.6  Numerical  Example  for  Approach  2 

Parameters  of  Example  1  (section  2.3.4)  are  used  for  this  numerical  example: 

Strut  lengths  an  =  3  in.  a34  =  3.5  in. 

Non-compliant  tie  lengths  a4i  =  4  in.  a23  =  2  in. 

Spring  1  free  length  &  spring  constant  L01  =  0.5  in.  ki  =  4  lbf/in. 

Spring  2  free  length  &  spring  constant  L02  =  1  in.  k2  =  2.5  lbf/in. 

For  this  example,  a  Maple  program  was  written  that  would  solve  the  determinant  of 

Sylvester  matrix  in  Equation  2-58.  The  result  was  a  polynomial  of  degree  32  with  respect 

to  the  variable  cos04.  The  solution  to  the  polynomial  for  this  numerical  example  resulted 

in  12  real  and  20  complex  roots.  Eight  of  the  real  roots  were  in  acceptable  range  of  -1  to 

+  1.  For  each  value  of  cos04  the  corresponding  value  of  cos0i  was  calculated  such  that 

Equations  2-37  and  2-55  are  simultaneously  satisfied.  These  values,  as  well  as  calculated 

values  for  Li  and  L2,  are  presented  in  Table  2-5. 

t  feasible  real  solutions. _ _ _ 

C4  (radian) _ ct  (radian) _ Li  (inches) _ L2  (inches) 

0.8144902900  0.2120649698  ~  -5.4853950884~  2.3332963547 

0.7083900124  0.1385596475  -5.3221641782  -2.9008756697 

0.9290901467  -0.9154296793  -1,7405998094  -1.4951507917 

0.8840460933  -0.9381755562  -1.5760033787  1.8699490332 

0.9046343262  -0.9312332602  1.6280054527  1.7088706403 

0.9434161072  -0.8970754537  1.8628443599  -1.3543814070 

0.6228159543  0.0544134519  5.1289299905  -3.2880318246 

0.8042767224  0.2077177765  5.4758767915  2.3937944296 

Complex  values  for  ci  that  correspond  to  each  complex  value  of  C4  were 

determined  and  are  presented  in  Table  2-6.  Although  these  solutions  are  not  realizable,  it 
is  important  to  show  that  since  the  C4,  ci  pairs  satisfy  the  two  equations,  no  extraneous 
solutions  were  introduced  in  the  variable  elimination  method. 


Table  2-5.Eighi 

Case _ 

1 


2  _ - 

3  - 

4  - 

5  - 

6  _ - 

7 
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Table  2-6.Twenty  complex  solutions. 


Case 

C4 

Cl 

1 

-1.2513573  -0.08533821 

1.0200425  +  0.0625503  1 

2 

-1.2513573  +  0.0853382  1 

1.0200425  -  0.0625503  1 

3 

-1.2371852  -0.05868661 

1.0164651  +0.03791361 

4 

-1.2371852  +  0.05868661 

1.0164651  -  0.03791361 

5 

-1.0662774  -0.02229991 

-1.0214122  +  0.01103171 

6 

-1.0662774  +  0.02229991 

-1.0214122  -0.01103171 

7 

-1.0406024  -0.0082765  1 

.1.0447401  +  0.01045821 

8 

-1.0406024  +  0.0082765  I 

-1.0447401  -0.01045821 

9 

-1.0322412  -0.0344383  I 

-1.0269253  +  0.0377413  1 

10 

-1.0322412  +  0.0344383  I 

-1.0269253  -0.0377413  1 

11 

-1.0087191  -0.00283591 

-1.1658773  +  0.03384621 

Table  2-6.  Continued. 


Case 

C4 

Cl 

12 

-1.0087191  +0.00283591 

-1.1658773  -0.03384621 

13 

0.3360373  -0.41274301 

-0.9379955  +  0.6334257  I 

14 

0.3360373  +  0.41274301 

-0.9379955  -  0.6334257  I 

15 

0.3571365  -0.10681791 

-0.9796431  +  0.3840155  I 

16 

0.3571365  +  0.10681791 

-0.9796431  -0.3840155  1 

17 

0.5105685  -0.0343981  I 

-1.0755053  -0.2990613  1 

18 

0.5105685  +  0.0343981  I 

-1.0755053  +  0.2990613  1 

19 

0.5418577-0.17049791 

-1.1122907  +  0.42951361 

20 

0.5418577  +  0.17049791 

-1.1122907  -0.42951361 

Lastly,  there  were  four  real  values  of  cos04  that  were  not  in  the  range  -1<  cos04  <  1 . 


These  values  are  listed  in  Table  2-7.  Corresponding  values  for  cos0i  could  not  be 
obtained  that  would  satisfy  both  Equations  2-37  and  2-48.  Thus  it  must  be  concluded  that 
four  extraneous  roots  were  introduced  in  this  second  solution  approach. 


Table  2-7.  Four  non-feasib 


e  real  solutions. 


Case 

C4  (radian) 

1 

-0.7460150  1018 

2 

-1 .0119084 

3 

-1.0061580 

4 

0.7460150  1018 
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2.4.7  Comparison  of  Results  for  Approaches  1  and  2 

The  numerical  example  from  Section  2.3.4  was  examined  using  this  second 
approach.  It  was  not  possible  to  symbolically  expand  the  determinant  of  the  coefficient 
matrix  in  Equation  2-58.  Rather  the  coefficients  of  the  polynomials  A,  B,  D,  and  Ci 
through  Cio  were  obtained  numerically  and  then  the  detenninant  was  expanded  to  obtain 
the  single  polynomial  equation  in  C4.  The  numerical  example  of  the  second  approach 
resulted  in  a  32nd  degree  polynomial  in  the  variable  C4.  Table  2-5  shows  the  resulting 
values  of  C4  that  solved  this  polynomial  as  well  as  the  corresponding  values  of  ci  such 
that  Equations  2-37  and  2-55  are  simultaneously  satisfied.  The  real  solutions  are 
identified  in  Table  2-5.  Values  of  Table  2-5  were  tested  in  the  differential  of  potential 
energy  to  identify  equilibrium  cases.  Analysis  showed  that  cases  3,  4,  5,  and  6  of  Table 
2-5  represent  equilibrium  configurations.  These  cases  are  identically  similar  to  cases 
presented  in  Table  2-1. 

Comparison  of  the  results  of  the  numerical  example  using  the  two  approaches 
outlined  in  this  chapter  showed  that  both  approaches  agree  on  the  correct  numerical 
answer.  Hence,  both  approaches  are  valid  in  finding  equilibrium  configuration  were 


closed  fonn  solution  for  this  structure  can  not  be  found. 


CHAPTER  3 

THREE  SPRING  PLANAR  TENSEGRITY 

3.1  Introduction 

This  chapter  considers  the  case  of  a  planar  tensegrity  system  with  one  non- 
compliant  member,  a4i,  three  compliant  members,  springs  Li,  L2,  L3,  and  two  struts,  a-,4 
andai2  (see  Figure  3-1).  In  Figure  3-1  the  two  struts  do  not  intersect  but  one  passes 
through  a  slit  cut  in  the  other  one.  The  objective  of  this  study  is  to  find  all  possible 
equilibrium  configurations  in  stable  condition  (minimum  potential  energy)  when  given 
the  lengths  of  the  struts  and  the  non-compliant  tie  together  with  the  free  length  and  spring 
constants  for  the  three  compliant  members. 

The  device  shown  in  Figure  3-1  is  a  two  degree  of  freedom  system.  Two 
parameters  must  be  specified,  in  addition  to  the  constant  mechanism  parameters,  in  order 
to  define  the  configuration  of  the  device.  These  two  parameters  will  be  referred  to  as  the 
descriptive  parameters  for  the  system.  One  obvious  set  of  descriptive  parameters  are  the 
angles  04  and  0i.  Considering  the  non-compliant  member  a4i  as  being  fixed  to  ground, 
specification  of  04  will  define  the  location  of  point  3.  Similarly,  specification  of  0i  will 
define  the  location  of  point  2. 

Two  approaches  to  solve  this  problem  are  presented  in  this  chapter.  Both  aim  to 
find  a  set  of  descriptive  parameters  that  minimize  the  potential  energy  in  the  system.  In 
the  first  approach,  the  lengths  of  the  compliant  members  Li  and  L2  are  chosen  as  the 
descriptive  parameters.  Derivatives  of  the  potential  energy  equation  are  obtained  with 
respect  to  Li  and  L2  and  values  for  the  descriptive  parameters  are  obtained  such  that  these 
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derivatives  are  zero,  corresponding  to  either  a  minimum  or  maximum  potential  energy 
state.  In  the  second  approach,  the  cosines  of  the  angles  04  and  0i  were  chosen  as  the 
descriptive  parameters.  The  cosines  of  the  angles  were  chosen  rather  than  the  angles 
themselves  in  the  hope  that  the  resulting  equations  would  be  simpler  in  that,  for  example, 
a  single  value  of  cos04  accounts  for  the  obvious  symmetry  in  solutions  that  will  occur 
with  respect  to  the  fixed  member  a4i. 


Figure  3-1.  Tensegrity  structure  with  two  struts,  a  non-compliant  (tie),  and  three  springs. 


3.2  First  Approach  Problem  Statement  (Descriptive  Parameters  Li,  L2) 


The  problem  statement  can  be  explicitly  written  as: 
Given:  ai2,  a34  lengths  of  struts, 

a4i  length  of  non-compliant  tie 


ki,  Loi  spring  constant  and  free  length  of  compliant  tie  between  points  4 
and  2 

k2,  L02  spring  constant  and  free  length  of  compliant  tie  between  points  3 
and  1 

k3,  L03  spring  constant  and  free  length  of  compliant  tie  between  points  2 
and  3 

Find:  Li 


L2 


length  of  spring  1  at  equilibrium  position, 
length  of  spring  2  at  equilibrium  position, 
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L3  corresponding  length  of  spring  3  at  equilibrium  position 


3.2.1  Development  of  Geometric  Constraint  Equation 

Figure  3-2  shows  the  nomenclature  that  is  used.  Li,  L2,  and  L3  are  the  extended 
lengths  of  the  compliant  ties  between  points  4  and  2,  points  3  and  l,  and  points  2  and  3. 


Figure  3-2.  Planar  tensegrity  structure. 

The  analysis  starts  by  considering  different  triangles  in  this  structure.  Figure  3-3  shows 
the  triangle  formed  by  side  a34,  L3,  and  L|. 


Figure  3-3.  Triangle  4-3-2. 
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A  cosine  law  for  this  triangle  can  be  written  as 


+  ^—  +  Lj  a34  cos04’=  — 
2  2  1  J4  4  2 


Solving  for  cos04'  yields 


t  2  _  t  2  _  2 

cos04’=  3  1  &34 


2Lj  a34 

Figure  3-4  shows  the  triangle  formed  by  side  a4i,  al2,  and  Li. 
triangle  can  be  written  as 


+  +  l  a41  cos04"=  3-1^- 

2  2  1  41  4  2 


Solving  for  cos04"  yields 


cos  04"  =  ai2'  L'~  a4i' 
2L,  a41 


2 


(3-1) 


(3-2) 

A  cosine  law  for  this 


(3-3) 


(3-4) 


Figure  3-4.  Triangle  4-1-2. 
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Figure  3-5  shows  the  triangle  formed  by  a4i,  a34,  and  L2.  A  cosine  law  for  this 
triangle  can  be  written  as 


34 


Ml 


+  a34  a41  cos04 


L, 


(3-5) 


3 


1-2 


1 


Figure  3-5.  Triangle  4-1-3. 
Solving  for  cos04  yields 


T  2  2  2 

0'-‘2  —  ^34  _  ^41 

4  =  — - — - — 

2^34  a41 


From  Figure  3-2  it  is  apparent  that 

04  +  04'  =  n  +  04" 

Equating  the  cosine  of  the  left  and  right  sides  of  Equation  3-7  yields 

COS  (04  +  04')  =  COS  (7T  +  04") 
and  expanding  this  Equation  3-8  yields 

cos04  cos04’  -  sin04  sin04'  =  -  cos04"  . 

Rearranging  (3-9)  gives 

cos04  cos04'  +  cos04"  =  sin04  sin04' . 


(3-6) 

(3-7) 

(3-8) 

(3-9) 

(3-10) 
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Squaring  both  sides  yields 

(cos04)2  (cos04')2  +  2  cos04  cos04'  cos04"  +  (cos04")2  =  (sin04)2  (sin04')2  (3-11) 

Substituting  for  (sin04)  and  (sin04')  in  terms  of  cos04  and  cos04'  gives 

(cos04)2  (cos04')2  +  2  COS04  cos04’  cos04"  +  (cos04")2  =  (1-cos204)  (1-cos204')  (3-12) 

Equations  3-2,  3-4,  and  3-6  are  substituted  into  Equation  3-12  to  yield  a  single 
equation  in  the  parameters  Li,  L2,  and  L3  which  can  be  written  as 


Gi  L34  +  (G2  L22  +  G3)  L32  +  (G4  l24  +  G5  L22  +  G6)  =  0 


(3-13) 


where 


Gi  -  a4]  , 

G2  =  G2aL12  +  G2b, 
G3  =  G3a  Li2  +  G3b, 


G4  -  G4a  Li  , 


G5  =  G5aL14  +  G5bL12  +  G5C, 


Gb  -  G6a  LE  +  G6b  . 


(3-14) 


and  where 


G2a  =  -1, 


2  2 
G2b  =  ai2‘  -  a4i  , 


2  2 

G3a  =  (a34  -  a4i  ), 


2  2  2  2  2  2 
G3b  =  a4i  (a4i  -  ai2  -  a34  )  -  ai2  a34  , 


G4a  =  1 , 


GSa=l, 


2  2  2  2  2  2 
Gst,  =  (-a  12“  -  a34“  -  a4i  ),  Gsc  =  a34  (a4i  -  ai2~), 


G6a  -  ai22  (a4i2  -  a342),  Gbb  -  ai2z  a34Z  (ai2z  +  a3/  -  a4iz) . 


2  „  2  2  ,  _  2  _2x 


(3-15) 


3.2.2  Development  of  Potential  Energy  Equations 

The  potential  energy  of  the  system  can  be  written  as 


U  — ^ki(Li  L01)  +  ^  k2  (L2  L02)  +  ^  k3  (L3  L03) 


(3-16) 
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At  equilibrium,  the  potential  energy  will  be  a  minimum.  This  condition  can  be 
determined  as  the  configuration  of  the  mechanism  whereby  the  derivative  of  the  potential 
energy  taken  with  respect  to  the  descriptive  parameters  Li  and  L2  both  equal  zero.  The 
geometric  constraint  equation.  Equation  3-13,  contains  three  unknown  terms,  Li,  L2,  and 
L3.  From  this  equation,  L3  can  be  considered  as  a  dependent  variable  of  Lj  and  L2.  The 
following  two  expressions  may  be  written: 


=ki(L1  L01)+k3(L3  L03)  3=0, 

oLj 

(3-17) 

dL 

=k2(L2  L02)+k3(L3  L03)  3  =0  . 

5L2 

(3-18) 

5L 

The  derivatives  8Lt,I8L\  and  SL3/SL2  can  be  determined  via  implicit  differentiation 
from  Equation  3-13  as 


3L3 

“Lj  [2Lj2  L22  +  G2a  L22  L32  +  G3a  L32  +  L24  +  G5b  L22  +  G6a] 

(3-19) 

3Lj 

L3[2G,  L32  +  G2a  Lj2  L22  +  G2b  L22  +  G3a  Lj2  +  G3b] 

0L3 

-  L2  [2Lj2  L22  +  G2a  Lj2  L32  +  G2b  L32  +  Lj4  +  G5b  Lj2  +  G5c] 

(3-20) 

8L2 

L3[2Gj  L32  +  G2a  Lj2  L22  +  G2b  L22  +  G3a  L,2  +  G3b] 

Substituting  (3-19)  into  (3-17)  and  rearranging  gives 

(Dj  L22+D2)  L33  +  (D3  L22+D4)  L32  +  (D5  L24+D6  L22+D7)  L3  +  (D8  L24+D9  L22+D10)  =  0 

(3-21) 


Dl  =  Dia  Li  , 

E>2  =  D2a  Li  +  D2b  , 
D3  =  D3a  Li  , 

D4  =  D4a  Li  , 

D5  =  E>5a  Li  , 


where 
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D6  =  Dga  Li3  +  D6b  Lj2  +  Dbc  Li  +  D6d  , 

D7  =  D7a  Li3  +  D7b  Li2  +  D7c  Li  +  D7d  , 

Ds  =  D8a  Li  , 

D9  =  Dga  L[3  +  D%  Li 
Dio  =  DioaLi  (3.22) 
and 

Dia  =  -  G2a  k3 , 

D2a  =  2  Gi  ki  —  G3a  k3,  D2b  =  -  2  Gi  ki  L01  , 

D3a  =  G2a  k3  L03 , 

D4a  =  G3a  k3  L03 , 

D5a  =  -  k3 , 

D6a  =  G2a  ki  -  2  k3 ,  D6b  =  -G2a  ki  L01  ,  D6c  =  G2b  ki  -  Gsb  k3 ,  D6d  =  -G2b  ki  L0i  , 

D7a =  G3a  ki  ,  D7b  =  -G3a  ki  L01  ,  D7c  =  G3b  ki  -  G6a  k3 ,  D7d  =  -G3b  ki  L01  , 

D8a  =  k3  L03 , 

Dga  =  2  k3  L03 ,  D%  =  G5b  k3  L03 

Dioa  =  G6a  k3  L03  (3-23) 

Substituting  (3-20)  into  (3-18)  and  rearranging  gives 
(El  L2  +  E2)  L33  +  (E3  L2)  L32  +  (E4  L23  +  E5  L22  +  E6  L2  +  E7)  L3  +  (E8  L23  +  E9  L2)  =  0 

(3-24) 


Ei  =  Eia  Lf  +  Eib  , 
E2  =  -  2  Gi  k2  L02  , 
E3  =  E3a  Lf  +  E3b  , 
E4  =  E4a  Lf  +  E4b  , 


where 
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E5  =  E5a  Lf  +  E5b  , 

E6  =  Eba  Li4  +  E6b  Li2  +  Egc , 

E7  =  E7a  Lr  +  E7b  , 

Es  =  Ega  Lj2  , 

E9  =  E9a  L[4  +  E9b  Lr  +  E9C  (3-25) 

and 

Ela  =  -G2a  k3  ,  Eib  =  -G2I5  k3  +  2  Gi  k2  , 

E3a  =  G2a  k3  L03  ,  E3b  =  G2b  k3  L03  , 

E4a =  G2a  k2  -  2  k3 ,  E4b  =  G2b  k2 , 

Esa  =  -G2a  k2  L02  ,  E5b  =  -G2b  k2  L02  , 

E6a  =  -  k3 ,  E6b  =  G3a  k2  -  G5b  k3 ,  Ebc  =  G3b  k2  -  Gsc  k3 , 

E7a  =  -G3a  k2  L02  ,  E7b  =  -G31,  k2  L02  , 

Esa  =  2  k3  L03  , 

E9a  =  k3  L03  ,  E9b  =  G5b  k3  L03  ,  Egc  =  Gsc  k3  L03  .  (3-26) 

3.2.3  Create  Solution  Matrix 

Equations  3-13,  3-21,  and  3-24  are  three  equations  in  the  three  unknowns  Li,  L2, 
and  L3.  Sylvester’s  method,  reference  1 1,  is  applied  in  order  to  obtain  sets  of  values  for 
these  parameters  that  simultaneously  satisfy  all  three  equations.  In  this  solution,  the 
parameter  Li  is  embedded  in  the  coefficients  of  the  three  equations  to  yield  three 
equations  in  the  apparent  unknowns  L2  and  L3.  Detennining  the  condition  that  these  new 
coefficients  (which  contain  Li)  must  satisfy  such  that  the  three  equations  can  have 
common  roots  for  L2  and  L3  will  yield  a  single  polynomial  in  Li. 
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2  2  2  2  2  2  3 

Equation  3-13  was  multiplied  by  L2,  L3,  L2L3,  L3  ,  L2  ,  L2L3  ,  L2  L3,  L2~L3  ,  L3  , 
L23,  L2L33,  L22L33,  L23L3,  L23L32,  and  L23L33.  Equation  3-21  was  multiplied  by  L3,  L2, 

L3",  L2L32,  L33,  L22,  L2L33,  L2  L3,  L22L32,  L23,  L34,  L23L3,  L23L3^,  and  L2L34.  Equation  3- 
24  was  multiplied  by  L2,  L3,  L2L3,  L32,  L22,  L2L32,  L22L3,  L22L32,  L33,  L23,  L2L33,  L22L33, 
L23L3,  L23L32,  L24,  L34,  L24L3,  L24L32,  L2L34,  and  L22L34.  This  resulted  in  a  set  of  52 
equations  that  can  be  written  in  matrix  for  as 

M).  =  0.  (3-27) 


The  vector  X  is  written  as 

X  =  [L2?L33,  L25L35,  L23L3?,  L2?L32,  L26L33,  L25L34,  L24L35,  L23L36,  L22L3?,  L2?L3, 

L26L32,  L25L33,  L24L34,  L23L35,  L2~L36,  L2L37,  L2?,  L26L3,  L25L3~,  L24L33,  L23L34, 
L2  L33,  L2L36,  L37,  L26,  L25L3,  L24L3^,  L23L33,  L2"L34,  L2L35,  L36,  L25,  L24L3, 

L23L32,  L2  L33,  L2L34,  L35,  L24,  L23L3,  L22L32,  L2L33,  L34,  L23,  L2L3,  L2L32,  L33, 


L22,  L2L3,  L32,  L2,  L3,  1]T  . 


(3-28) 


The  coefficient  matrix  M  is  a  52^52  matrix.  To  show  this  52x52  coefficient 
matrix,  it  is  subdivided  into  the  four  26x26  submatrices  Mu,  Mi2,  M2i,  M22  as 


VI 


Mil  Ml2 
M21  M22 


(3-29) 


The  four  sub  matrices  are: 
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'00000000000  G40G20G,  ooo  og5og3o  og6 

0  0  0  0  0  0  D5  0  D,  0  0  D8  0  D3  0  0  0  D6  o  d2  d9  0  d4  0  d7  d10 

0  0  0  0  0  0  0  0  0  o  o  o  e4  o  e,  o  e8  e5  e3  e2  o  e6  o  e9  e7  0 

D5  0  D,  0  0  0  D8  0  D3  0  0  0  0  D6  0  D2  0  D9  0  D4  0  0  D7  0  D]0  0 

0  0  0  0  0  0  E4  0  E,  0  o  e8  e5  e3  e2  0  0  e6  0  o  e9  e7  0  0  0  0 

0  0  0  0  0  0  0  E4  0  E,  0  o  e8  e5  e3  e2  0  0  e6  0  o  e9  e7  0  0  0 

E4  0  E,  0  0  0  E8  E5  E3  E2  0  0  0  E6  0  0  0  E9  E7  0  0  0  0  0  0  0 

00000  g4  0  g2  0  G,  0000  0  0  G50G30  ooog6o  0 

000000  G40G20G,  000  0  00  G50G3  0  0  0  0  G60 

0  g2  0G,  00000000  g5og3  000  0  0  0  G6  0  0  0  0 

0  D,  0  0  0  D8  0  D3  0  0  0  0  D6  0  D2  0  D9  0  D4  0  0  D7  0  D10  0  0 

D8  0  D3  0  0  0  0  0  D6  0  D2  0  0  D9  0  D4  0  0  0  D7  0  0  Dl0  0  0  0 

_  0  d6  o  d2  ooo  d9  o  d4  0  0  0  0  d7  0  0  0  D]0  0  0  0  0  0  0  0 

12  “  0  E4  0  E,  0  0  o  e8  e5  e3  e2  0  0  0  E6  0  0  0  e9  E  0  0  0  0  0  0 

0  E,  0  0  o  e8  e5  e3  e2  0  0  0  E6  0  0  o  e9  e7  0  0  0  0  0  0  0  0 

E8  E5  E3  E2  0000  E6  0000  E9E7  000  0  0  OOO  0  0  0 

E5  E3  E2  0000  E6  0000  E9E7  0  000  0  0  OOO  0  0  0 

0E6  00000  E9E7  0000000000  0  OOO  0  0  0 

G4  0  G2  OG,  00000000  G5OG3  OOO  0  0  0  G6  0  0  0 

G2  0  G,  00000000  G5OG3  0  000  0  0  G6  0  0  0  0  0 

0000000  G5OG3  00000000  G60  OOO  0  0  0 

000000  G50G3  00000000  G6  0  0  OOO  0  0  0 

G5OG3  0000000000  G6  00000  0  OOO  0  0  0 

00D60D2  000  D90D4  000  0  D7  00  ODoO  0  0  0  0  0 

D3  0  0  0  0  0  D6  0  D2  0  0  D9  0  D4  0  0  0  D7  0  0  D10  0  0  0  0  0 

o  d9  od4  00000  D7  0000  D10  000  0  0  OOO  0  0  0 

D6  0  D2  0  0  0  D9  0  D4  0  0  0  0  D7  0  0  0  D10  0  0  0  0  0  0  0  0 

d9  o  d4  0  0  0  0  0  D7  0  0  0  0  D10  0  0  0  0  00  0  0  0  0  0  0 

0E8E5E3E2  00  0  0  E6  000  oe9e7o  0  00  0  0  0  0  0  0 

E3  E2  0000  E6  0000  E9E7  OOO  0  0  OOO  0  0  0  0  0 

00  E6  00000  E9E7  OOO  OOO  0  0  OOO  0  0  0  0  0 

0E9E7  0000000000  OOO  0  0  OOO  0  0  0  0  0 

E6  00000  E9E7  0  0000  00  0  0  0  00  0  0  0  0  0  0 

E9E7  0000000000000000  00000000 

00000000  G50G3  00  OOO  0  0  0  G6  0  0  0  0  0  0 

00000  G50G3  00000  000  G60  OOO  0  0  0  0  0 

OG5OG3  000000000  0  G6  0  0  0  OOO  0  0  0  0  0 

00000000  G6  0000  00  0  0  0  00  0  0  0  0  0  0 

0G3  0000000000  G6  00000  00000000 

22  "0000000  G6  00000  OOO  0  0  00  0  0  0  0  0  0 

0G6  0000000000000000  00000000 

o  d2  ooo  d9  o  d4  0  0  0  0  d7  ooo  D10  0  OOO  0  0  0  0  0 

00D90D4  00000  D7  00  ood10o  0  00  0  0  0  0  0  0 

0D4  00000  D7  0000  D10  OOO  0  0  OOO  0  0  0  0  0 

0D7  00000  D10  00000  OOO  0  0  OOO  0  0  0  0  0 

OOO  D7  000  0  0  D10  OOO  OOO  0  0  OOO  0  0  0  0  0 

00000  E9E7  000000  OOO  0  0  00  0  0  0  0  0  0 

OOO  E6  00000  E9E7  00  OOO  0  0  OOO  0  0  0  0  0 

E7  0000000000000000000000000 
00000000000000000000000000 
00  E9E7  000000000  OOO  0  0  OOO  0  0  0  0  0 

00000000000000000000000000 


(3-30) 
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In  order  for  a  solution  for  the  vector  X  to  exist,  it  is  necessary  for  the  52  resulting 
equations  to  be  linearly  dependent.  This  will  occur  if  the  determinant  of  M  equals  zero. 

It  was  not  possible  to  symbolically  expand  the  determinant  of  matrix  M.  A 
numerical  case  was  analyzed  and  a  polynomial  of  degree  158  in  the  variable  Li  was 
obtained.  It  was  not  possible  to  solve  this  high  degree  polynomial  for  the  values  of  Li, 
although  several  commercial  and  in-house  written  algorithms  were  attempted.  Because 
of  this,  a  different  method  was  attempted  to  solve  the  set  of  equations  3-13,  3-21,  and 
3-24. 


3.2.4  Solution  of  Three  Simultaneous  Equations  in  Three  Unknowns  - 
Continuation  Method 

The  continuation  method  (Morgan  [18],  Mora  [19],  Cullen  [20],  Strang  [21],  Garcia 
and  Li  [22],  Morgan  [23-24],  Wampler  et  al.,  [25])  which  is  a  numerical  technique  to 
solve  a  set  of  equations  in  multiple  variables  are  used.  This  is  as  opposed  to  Sylvester’s 
method  which  would  lead  to  a  symbolic  solution  of  the  problem. 

A  concise  description  of  the  continuity  method  is  presented  by  Tsai,  1999. 

Suppose  one  wishes  to  solve  the  set  of  equations  F(x)  which  are  defined  by 


F(x): 


fj(xi,x2,---,xn)  =  0 
f2  (x1,x2,---,xn)  =  0 


lf„(Xl>X2,'--,Xn)  =  0 


(3-31) 


F(x)  is  called  the  target  system. 

The  continuation  method  begins  by  first  estimating  the  total  number  of  possible 
solution  sets  (sets  of  values  for  Li,  L2,  and  L3  for  our  case)  that  satisfy  the  given 
equations.  For  example,  Bezouf  s  theorem  states  that  a  polynomial  of  total  degree  n  has 
at  most  n  isolated  solutions  in  the  complex  Euclidean  space.  Including  solutions  at 
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infinity,  the  Bezout  number  of  a  polynomial  system  is  equal  to  the  total  degree  of  the 
system. 

Next,  an  initial  system,  G(x)  =0,  is  obtained,  whose  solution  will  be  of  the  same 
degree  as  that  of  F(x),  but  whose  solution  set  is  known  in  closed  form.  In  other  words, 
G(x)  maintains  the  same  polynomial  structure  as  F(x). 

Finally,  a  homotopy  function  H(x,  t)  is  prepared  such  as 

H(x,  t)  =  y  (1-t)  G(x)  + 1  F(x)  (3-32) 

where  y  is  a  random  complex  constant.  When  t=0,  the  homotopy  function  equals  the 
initial  system,  G(x).  When  t=l,  the  homotopy  function  equals  the  target  system,  F(x). 
Recall  that  the  solutions  to  G(x)  are  known.  As  the  parameter  t  is  increased  in  small  steps 
from  0  to  1,  the  solutions  of  H(x,  t)  can  be  tracked  (referred  to  as  path  tracking)  and  when 
t  =1,  these  solutions  will  be  the  solutions  to  the  original  target  system.  If  the  degree  of 
the  solution  set  was  overestimated,  some  of  the  solutions  will  track  to  infinity  and  these 
can  easily  be  discarded. 

3.2.5  Numerical  Example 

The  following  infonnation  is  given: 
strut  lengths: 

an  =  14  in.  a34  =  12  in. 
non-compliant  tie  lengths: 

a4i  =  10  in. 

spring  1  free  length  &  spring  constant: 

Loi  =  8  in.  ki  =  1  lbf/in. 
spring  2  free  length  &  spring  constant: 


L02  =  2  in. 


k2  =  2.687  lbf/in. 
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spring  3  free  length  &  spring  constant: 

L03  =  2.5  in.  k3  =  3.465  lbf/in. 

Based  on  these  values,  the  coefficients  in  equations  3-13,  3-21,  and  3-24  were 
evaluated  to  yield  the  three  equations 

100  L34  +  [(-Li2  +  96)  L22  +  44  Lj2  -  52224]  L32  +  L\  L24 
+  (Li4  -  440  Li2  -13824)  L22  -  8624  Li2  +  6773760  =  0  (3-33) 

(2.5  L|  L22  +  90  Li  -  1600)  L33  +  (-8.663  Li  L22  +  381.165  Lj)  L32  +  [-2.5  L,  L24 
+  (-6  Li3  +  8  Li2  +  1 196  Li  -  768)  L22  +  44  Lj3  -  352  Ll2  -  30664  Li  +  417792]  L3 

+  8.663  L,  L24  +  (17.326  Lj3  -  381 1.651  U)  L22  -  74708.354  Li  =  0  (3-34) 

[(2.5  Lf  +  160)  L2  -  1074.637]  L33  +  (831.633  -  8.663  L!2)  L2  L32 
+  [(-7  Li2  +  192)  L23  +  (-515.826  +  5.373  Li2)  L22  +  (-2.5  Lj4  +  1188  Lx2  -  69888)  L2 
-  236.420  Li2  +  280609.161]  L3  +  17.326  Lj2  L23 
+  (-119755.135  +  8.663  Li4- 3811.651  Lj2)  L2  =  0  .  (3-35) 

The  continuation  method  was  run  on  this  set  of  three  equations  in  three  unknowns 
to  obtain  all  solution  sets  for  the  three  spring  lengths,  Li,  L2,  and  L3,  for  the  particular 
numerical  example.  The  software  PHCpack  (Verschelde  [26])  was  used  to  implement  the 
method. 

The  PHCpack  software  estimated  the  number  of  possible  solutions  to  be  136. 


Seven  real  solutions  were  obtained  and  these  are  listed  in  Table  3-1. 
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Table  3-1:  Seven  real  solutions  for  three-spring  planar  tensegrity  system(units  are  inches) 
Case  Li  L2  L3 


1 

13.000 

8.000 

7.017 

2 

-11.376 

-10.371 

-5.333 

3 

-7.585 

9.097 

10.106 

4 

-11.029 

12.557 

-3.044 

5 

13.969 

-5.800 

9.164 

6 

14.248 

-9.373 

-4.774 

7 

13.181 

11.599 

-2.488 

An  equilibrium  analysis  was  conducted  for  the  seven  cases  and  only  the  first  case 
was  indeed  in  equilibrium.  Case  1  is  shown  in  Figure  3-6. 


Figure  3-6.  Case  1,  equilibrium  solution. 

3.3  Second  Approach  Problem  Statement  (calculate  cos04,  cosOi) 

In  the  second  approach,  the  cosines  of  the  angles  04  and  9i  were  chosen  as  the 
descriptive  parameters.  The  cosines  of  the  angles  were  chosen  rather  than  the  angles 
themselves  in  the  hope  that  the  resulting  equations  would  be  simpler  in  that,  for  example, 
a  single  value  of  cos04  accounts  for  the  obvious  symmetry  in  solutions  that  will  occur 
with  respect  to  the  fixed  member  azn. 

The  problem  statement  is  presented  as  follows: 
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Given:  Strut  lengths  ai2,  a34,  tie  length  a4i,  and  the  spring  constants  and  free 
lengths  ki,  k2,  k3,  L0i,  L02,  L03. 


Find:  cos04  and  cos0i  when  the  system  is  in  equilibrium. 

The  solution  approach  is  as  follows: 

1 .  Obtain  expressions  for  Li,  L2,  and  L3  in  terms  of  cos04  and  cos0i 

2.  Write  the  potential  energy  equation 

3.  Determine  values  of  cos04  and  cos0i  such  that  dU/dcos04  =  dU/d  cos0i=  0 
3.3.1  Development  of  Geometric  Equations 

Figure  3-1  shows  the  nomenclature  used.  Springs  Li,  L2,  and  L3  are  the  extended 
lengths  of  the  compliant  ties  between  points  4  and  2,  2  and  3,  and  1  and  3  respectively. 
The  unloaded  lengths  of  springs  are  given  by  L0i,  L02,  and  Lo3.  A  cosine  law  for  the 
quadrilateral  1 -2-3-4  can  be  written  as: 


where, 


and  where 


2 

z41  =  — ai2  (X4Si  +  Y4Cj  )  +  Z4  4  — 


X4  -  a34  S4 
Y4  =  -  (a4i  4-  a34  C4) 


2  2 

Z&34  .  3.41 

4  -  +  +  a34a41C4  ' 


(3-36) 


(3-37) 


(3-38) 

(3-39) 

(3-40) 


Substituting  (3-37)  into  (3-36)  and  rearranging  gives 
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3.J2X4S1  “I-  Z4 


*12 


(3-41) 


Substituting  (3-38),  (3-39),  and  (3-40)  into  (3-41),  then  squaring  it,  and  substituting  for 

2  2  22 

S4  =  1-  C4  and  Si  =l-cf  and  multiplying  the  entire  equation  by  4  yields 

L34  +  A  L32  +  B  =  0  (3-42) 

where 


A  —  Ai  C4  +  A2  Ci  +  A3  Ci  C4  +  A4  ,  (3-43) 

B  =  Bi  cE  +  B2  cr  C4  +  B3  Ci  +  B4  Ci  C4  +  B5  Ci  04“  +  B6  C4”  +  B7  C4  +  Bg  (3-44) 


and  where 


Ai  =  -  4  a34  a4i 
A2  =  -  4  ai2  &4i 
A3  =  -  4  ai2  a34 
A4  =  -  2  (an2  +  a342  +  a4i2) 

Bi  =  4  ai22  (a342  +  a4i2) 

B2  =  8  ai2~  a34  a4i 

B3  =  4  ai2  a4i  (a4i“  +  ai22  +  a342) 

B4  =  4  ai2  a34  (a34_  +  ai2  +  3  a4i  ) 

B5  =  8  ai2  a4i  a34“ 

B6 =  4  a34*’  (an2  +  a4i") 

B7  =  4  a34  a4i  (a4i“  +  a34~  +  ai22) 

Bg  =  (a4i2  +  a342  -  2  ai2  a34  +  ai22)  (a4f  +  a342  +  2  ai2  a34  +  ai22)  .  (3-45) 

Equation  3-42  expresses  L3  as  a  function  of  C4  and  ci.  A  cosine  law  for  triangle  4- 


1  -2  may  be  written  as 
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&19  ^41  -L/1 

-^  +  ^y-  +  apa41c1  =  — 
2  2  2 


(3-46) 


A  cosine  law  for  the  triangle  3-4-1  may  be  written  as 


2  2  j  2 

^34  .  ^41  .  ^-'7 

2  2  34a41C4  =-J- 


(3-47) 


Equations  3-42,  3-46,  and  3-47  define  the  three  spring  lengths  in  terms  of  the  variables 
cos04  and  cosOi. 


3.3.2  Development  of  Potential  Energy  Equations 

The  total  potential  energy  stored  in  all  three  springs  is  given  by, 

U  =  >/2  k!  (Li-Loi)2  +  Ei  k2  (L2-L02)2  +  V2  k3  (L3-L03)2  .  (3-48) 

Differentiating  the  potential  energy  with  respect  to  C4  and  ci  and  then  evaluating  values 
for  C4  and  ci  that  cause  the  derivative  of  the  potential  energy  to  equal  zero,  will  identify 
configurations  of  either  minimum  or  maximum  potential  energy.  These  derivatives  may 
be  written  as, 

dU/dc4  =  ki  (Li  -  L01)  dLi/dc4  +  k2  (L2  -  L02)  dL2/dc4  +  k3  (L3  -  L03)  dL3/dc4 

(3-49) 

dU/dci  =  ki  (Li  -  L01)  dLi/dci  +  k2  (L2  -  L02)  dL2/dci  +  k3  (L3  -  L03)  dL3/dci 

(3-50) 

Since  from  3-46,  Li  is  not  a  function  of  C4,  dLi/dc4  =  0.  Similarly,  from  3-47,  L2  is  not  a 
function  of  ci  and  thus  dL2/dci  =  0.  Equations  3-49  and  3-50  reduce  to 

dU/dc4  =  k2  (L2  -  L02)  dL2/dc4  +  k3  (L3  -  L03)  dL3/dc4  ,  (3-51) 

dU/dci  =  ki  (Li  -  L01)  dLi/dci  +  k3  (L3  -  L03)  dL3/dci  .  (3-52) 

The  term  dL2/dc4  is  evaluated  from  3-47  as 
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^2  _  a34a41 

dc4  L2 

and  the  term  dLi/dci  is  evaluated  from  3-46  as 

dLt  _  a12a41 
dc!  Lj 

Implicit  differentiation  of  3-42  for  L3  with  respect  to  C4  and  ci  yields 

dL-3  (A[  A3Cj)L3  +  ByCj  ®4^1  2^B5C  1C4  +  2B6C4  +  By 

dc4  2L3(2L3  +  A,c4  +  A2CJ  +  A3c1c4  +  A4) 


(3-53) 


(3-54) 


(3-55) 


dL3  (At  +  A3c4)L3  +  2B,Cj  +  2B-,c^c4  +  B3  +  B4c4  +  B3c4” 
dCj  2L3(2L3  +  AjC4  +  AtC^  +  A3C[C4  +  A4) 


(3-56) 


Substituting  3-53  and  3-55  into  3-51  and  3-54  and  3-56  into  3-52  and  equating  to  zero 
yields 

(Mi  L33  +  M2  L32  +  M3  L3  +  M4)  L2  +  M5  L33  +  M6  L3  =  0  (3-57) 

(N  ,  L33  +  N2  L32  +  N3  L3  +  N4)  L!  +  N5  L33  +  N6  L3  =  0  (3-58) 

Where  the  coefficients  Mj  and  N,,  i=1..6,  are  functions  of  c4  and  Ci  as 
Mi  =  -k3  (A3  ci  +  Ai)  +  4  k2  a34  a4i 
M2  =  k3  L03  (A3  ci  +  Ai) 

M3  =  -k3  B2  ci2  +  [2  (k2  a34  a4iA3  -  k3  B5)  c4  +  2  k2  a34  a4i  A2  -  k3  B4)  Ci 
+  2  (k2  a34  a4i  Ai  -  k3  B6)  c4  -  k3  By  +  2  k2  a34  a4)  A4 
M4  =  k3  L03  B2  ci2  +  k3  L03  B4  Ci  +  2  k3  L03  B5  Ci  c4  +  2  k3  L03  B6  c4  +  k3  L03  B7 
M5  =  -4  k2  a34  a4i  Lo2 

M6  =  -2  k2  (a34  a4iLo2  A2  Ci  +  a34  a4i  Lo2  A3  Ci  c4  +  a34  a4i  Lo2  Ai  c4 

+  a34  a^ii  Lo2  A4)  (3-59) 

Ni  =  -  k3  (A3  c4  +  A2)  +  4  ki  ai2  a4i 
N2  =  k3  L03  (A3  c4  +  A2) 
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N3  =  2  ki  ai2  &4i  (Ai  C4  +  A2  Ci  +  A3  Ci  C4  +  A4)  —  k3  (B4  C4  +  B5  C4”  +  2  Bi  Ci 


+  2  B2  Ci  C4  +  B3) 

N4  =  k3  L03  (B4  C4  +  B5  C4”  +  2  Bi  Ci  +  B2  Ci  C4  +  B3) 

N5  =  -  4  ki  ai2  a4i  L01 

N6  =  -2  ki  ai2  a4i  L0i  (A4  +  Ai  c4  +  A2  Ci  +  A3  ci  c4)  .  (3-60) 

Equations  3-57  and  3-58  may  be  written  as 

(Mi  L33  +  M2  L32  +  M3  L3  +  M4)  L2  =  -  M5  L33  -  M6  L3 ,  (3-61) 

(Ni  L33  +  N2  L32  +  N3  L3  +  N4)  Li  =  -  N5  L33  -  N6  L3 .  (3-62) 

Squaring  both  sides  of  both  equations  gives 

(M,  L33  +  M2  L32  +  M3  L3  +  M4)2  L22  =  (Ms  L33  +  M6  L3)2  ,  (3-63) 

(Ni  L33  +  N2  L32  +  N3  L3  +  N4)2  Li2  =  (Ns  L33  +  N6  L3)2  .  (3-64) 


2  2 

Using  3-46  and  3-47  to  substitute  for  Li  and  L2~  will  yield  two  equations  in  the 
parameters  Ci,  C4,  and  L3.  These  two  equations  can  be  arranged  as 

(Pi  Ci2  +  p2  Ci  +  p3)  C43  +  (p4  Ci3  +  p5  Ci2  +  p6  Ci  +  p7)  C42 
+  (p8  Cl4  +  p9  Cl3  +  P10  Cl2  +  P11  Cl  +  P12)  C4 
+  (p  1 3  Cl4  +  P14  Cl3  +  p  1 5  Cl2  +  Pi6  Cl  +  P17)  =  0  (3-65) 

(qi  Ci  +  q2)  c44  +  (q3  Ci2  +  q4  Ci  +  q5)  c43 
+  (q6  Ci3+  q7  cr  +  q8  Ci  +  q9)  C4~  +  (qio  Ci3  +  qn  cf  +  qi2  ci  +  qn)  C4 

+  (qi4  Ci3  +  qis  cf  +  qi6  Ci  +  qn)  =  0  (3-66) 

where 

Pl  =PlaL32  +  Plb  L3  +  pic 
P2  =  P2a  L32  +  p2b  L3  +  p2c 
P3  =  P3a  L32  +  p3b  L3  +  p3c 
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P4  —  P4a  L3  +  P4b  L3  +  P4c 

P5  =  P5a  L34  +  p5b  L33  +  P5c  ^  +  P5d  L3  +  p5e 

P6  P6a  L3  Pbb  L3  Pbc  L3  Pbd  L3  +  P6e 

p7  =  P7a  L34  +  p7b  L33  +  p7c  L32  +  p7d  L3  +  P7e 

P8  =  P8a  L32  +  P8b  L3  +  P8c 

P9  —  P9a  L3  +  p%  L3  +  P9c  L3  P9d  L3  +  P9e 

PlO  —  PlOa  L3  +  PlOb  L3  +  PlOc  L3  +  PlOd  L3  +  PlOe  L3  +  PlOf  L3  +  PlOg 

Pi  1  =  Pi  la  L36  +  Pi  lb  L35  +  Pile  L34  +  Pi  Id  L33  +  Pile  L32  +  pi  if  L3  +  pi  ig 

P12  =  Pl2a  L36  +  Pi 2b  L35  +  Pl2c  I-34  +  Pl2d  I-33  +  Pl2e  1-3“  +  Pl2f  L3  +  Pl2g 

2 

P 1 3  =  Pl3a  1-3“  +  Pi  3b  1-3  +  Pl3c 

4  3  2 

Pl4  —  Pl4a  L3  +  Pl4b  L3  +  Pl4c  L3  +  Pl4d  L3  +  Pl4e 

P 15  =  Pl5a  L36  +  Pi 5b  I-35  +  Pi 5c  I-34  +  Pl5d  I-33  +  Pl5e  1-3“  +  Pl5f  1-3  +  Pl5g 

Pl6  Pl6a  L3  Pl6b  L3  Pl6c  L3  Pl6d  L3  Pl6e  L3 ”  Pl6f  L3  Pl6g 

P 1 7  =  Pl7a  L36  +  Pi 7b  I-35  +  Pl7c  L34  +  Pnd  L33  +  pi7e  hf  +  Pl7f  L3  +  Pi7g 

qi  =  qia  L3”  +  qib  L3  +  qic 

q2  =  q2a  1-3“  +  q2b  L3  +  q2c 

q3  =  q3a  1-3“  +  q3b  L3  +  q3C 

q4  —  q4a  L3  +  q4b  L3  +  q4C  L3  +  q4<i  L3  +  q4e 

qs  =  qsa  L34  +  qsb  L33  +  qsc  +  qsd  L3  +  qse 

q6  q6a  L3  +  q(,h  L3  +  q6c 

q7  =  q7a  L34  +  q7b  L33  +  q7c  L32  +  q7d  L3  +  q7e 

qs  =  qsa  L36  +  qsb  L35  +  q8c  L34  +  q8d  L33  +  q8e  L32  +  q8f  L3  +  q8g 

q9  =  q9a  L36  +  q9b  L35  +  q9C  L34  +  q9d  L33  +  q9e  L32  +  q9f  L3  +  q9g 


(3-67) 
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qio  =  qioa  L3"  +  qiob  L3  +  qioc 

q  1 1  —  qua  L3  +  qnb  L3  +  qiic  L3  +  qiid  L3  +  qiie 

qi2  =  qi2a  L36  +  qi2b  L35  +  qi2c  L34  +  qi2d  L33  +  qi2e  +  qi2f  L3  +  qi2g 

qi3  =  qi3a  L36  +  qi3b  L35  +  qnc  L34  +  qi3d  L33  +  qi3e  L3”  +  qi3f  L3  +  qi3g 

qi4  =  qi4a  L3"  +  qi4b  L3  +  qi4c 

q  1 5  —  qi5a  L3  +  qisb  L3  +  qi5c  L3  +  qisd  L3  +  qise 

qi6  qi6a  L3  qi6b  L3  qi6C  L3  qi6d  L3  -t-  qi6e  L3”  qi6f  L3  qi6g 

qi7  =  qi7a  L36  +  qi7b  L35  +  qi7c  L34  +  qm  L33  +  qi7e  +  qi7f  L3  +  qng  •  (3-68) 

The  coefficients  pia  through  qi7g  are  functions  of  the  given  constant  parameters.  These 
coefficients  are  defined  in  the  Appendix. 

The  coefficients  pia  through  qi7g  have  been  obtained  symbolically.  For  example, 
the  terms  pia  through  pic  are  written  as 

Pia  =  8  k22  a343  a4i3  A32  -  16  k2  a342  a4]2  A3  k3  B5  +  8  k32  B52  a34  a4i 
pib  =  -16  k3_  B52  L03  a34  a4i  +  16  k2  a34“  a4f  A3  k3  L03  B5 
Pic  =  8  k3“  L03"  B5"  a34  a4i 

The  remaining  terms  are  not  listed  here  due  to  their  complexity. 

Similarly,  equation  3-42  is  an  equation  in  the  same  parameters,  Ci,  c4,  and  L3.  This 
equation  can  be  factored  as 

(ri  Ci  +  r2)  c42  +  (r3  Ci2  +  r4  Ci  +  r5)  c4  +  (r6  Cj2  +  r7  Ci  +  r8)  =  0  (3-69) 

rj  =  8  ai2  a4i  a34" 
r2  =  4  a34"  (ai22  +  a4i2) 
r3  =  8  a34  a4i  ai22 


where 
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r4  =  (-4  ai2  a34)  L32  +  4  a32  a34  (a342  +  ai22  +  3  a4i2) 


rs  -  (-4  a34  a4i)  L3“  +  4  a34  a4i  (a34  +  ai2 


r6  =  4  an2  (a342  +  a4i2) 


r7  =  (-4  ai2  a4i)  L3  +  4  a12  a4i  (a4f  +  a34‘ 


r8  =  L34  -  2  (a342  +  a122  +  a412)  L32  (3  -70) 

3.3.3  Solution  of  Equation  Set  for  c4,  Ci,  and  L3 

Equations  3-65,  3-66,  and  3-69  are  factored  such  that  the  parameter  L3  is  embedded 
in  the  coefficients  pi  through  rs.  Sylvester’s  method  was  used  in  an  attempt  to  solve 
these  equations  for  all  possible  sets  of  values  of  L3,  ci,  and  c4.  Equation  3-65  was 
multiplied  by  ci,  Ci  ,  c4,  c4  ,  c4  ,  ci  ,  Ci2e4,  CiC4  ,  and  ci  c4“  to  obtain  10  equations 

2  3  2  2  3 

(including  itself).  Equation  3-66  was  multiplied  by  ci,  cf,  ci  ,  c4,  c4  ,  cic4,  cfc4,  ci  c4, 

2  2  2* 

CiC4“  and  ci  c4  to  obtain  1 1  equations  (including  itself).  Equation  3-69  was  multiplied  by 

234  234  2  3  4  2223242  3 

Cl?  Cl  ,  Cl  ,  Cl  ,  C4,  C4  5  C4  9  C4  ,  C1C4,  Ci  C4,  Ci  C4,  Ci  C4?  C1C4  ,  Ci  C4  ,  Ci  C4  ?  Ci  C4  ,  C1C4  9 


Ci2c43,  Ci3Ci4,  CiC44  and  Ci2e44  to  obtain  22  equations  including  itself.  This  resulted  in  a  set 
of  43  “homogeneous”  equations  in  43  unknowns. 

The  condition  for  a  solution  to  this  set  of  equations  was  that  they  are  linearly 
dependent,  i.e.  the  determinant  of  the  43x43  coefficient  matrix  must  equal  zero.  This 
would  yield  a  polynomial  in  the  single  variable  L3.  Due  to  the  complexity  of  the 
problem,  it  was  not  possible  to  expand  this  determinant  symbolically.  Similarly,  it  was 
not  possible  to  root  the  polynomial  that  resulted  from  a  particular  numerical  example. 
Because  of  this,  the  Continuation  Method  was  used  again  to  determine  solutions  for  the 
same  numerical  example  presented  in  Section  3.2.5. 
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Substituting  the  given  numerical  values  into  equations  3-65,  3-66,  and  3-69  are 
yield  three  equations  for  this  numerical  example.  For  this  case,  the  coefficients  pia 
through  qi7g  and  ri  through  r8  are  evaluated  and  presented  in  Appendix  b. 

The  continuation  method  was  run  on  this  set  of  three  equations  in  three  unknowns 
to  obtain  all  solution  sets  for  the  three  variables,  L3,  Ci,  and  C4,  for  the  particular 
numerical  example.  The  software  PHCpack  (Verschelde  [26])  was  again  used  to 
implement  the  method. 

The  PFICpack  software  estimated  the  number  of  possible  solutions  to  be  136. 
Seven  real  solutions  were  obtained  and  these  are  listed  in  Table  3-2. 


Table  3-2.  Real  value  answers  using  cosine  of  Oi,  64  in  radians  and  L3  in  inches. 


Case 

Cosine  0i 

Cosine  04 

l3 

1 

-0.453571523 

-0.749999925 

7.01658755 

2 

-0.594965126 

-0.568509675 

-5.33336481 

3 

-0.851695156 

-0.671871253 

10.1059910 

4 

-0.622722694 

-0.359683532 

-3.04401046 

5 

-0.360201938 

-0.876513972 

9.16391281 

6 

-0.332122619 

-0.650641695 

-4.77378239 

7 

-.0436671807 

-0.456122910 

-2.48758125 

The  lengths  of  the  springs  Li  and  L2  were  calculated  for  each  of  the  seven  cases  and 


it  was  determined  that  the  seven  real  solutions  here  correspond  directly  with  the  seven 
real  solutions  that  were  detennined  from  the  first  problem  formulation  and  listed  in  Table 
3-1.  As  before,  case  1  was  the  only  case  that  was  indeed  in  equilibrium  and  the 
mechanism  is  shown  in  this  configuration  in  Figure  3-6. 

3.4  Conclusion 

Two  approaches  were  presented  to  solve  the  three-spring  planar  tensegrity 
equilibrium  problem.  Both  approaches  resulted  in  high  degree  polynomials.  These 
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polynomials  could  not  be  derived  symbolically  and  were  in  fact  difficult  to  derive 
numerically. 

A  continuation  approach  was  presented  to  solve  the  set  of  three  equations  for  each 
of  the  two  solution  approaches.  A  force  balance  analysis  was  conducted  to  identify  all 
realizable  solutions;  in  this  case  one. 

It  is  apparent  that  further  analysis  of  Sylvester’s  variable  elimination  procedure 
may  be  necessary  to  reduce  the  number  of  simultaneous  equations,  and  thereby  the  size 
of  the  univariate  polynomial  that  results.  As  a  result  of  the  work  in  this  chapter  it  can  be 
said  that  the  degree  of  the  solution  has  been  bounded,  although  further  work  may  yield  a 
simpler  result. 


CHAPTER  4 

FOUR  SPRING  PLANAR  TENSEGRITY 


4.1  Introduction 

In  this  chapter  planar  tensegrity  mechanisms  with  two  struts  and  four  elastic  ties 
(see  Figure  4-1)  are  analyzed  to  determine  all  possible  equilibrium  configurations  for  the 
device  when  no  external  forces  or  moments  are  applied.  The  stable  equilibrium  position 
is  determined  by  identifying  the  configurations  at  which  the  potential  energy  stored  in  the 
four  springs  is  a  minimum. 


Figure  4-1.  Two  strut,  four  spring  planar  tensegrity  device 

Four  non-linear  equations  are  obtained  in  the  variables  Li,  L2,  L3,  and  L4  which 
correspond  to  the  equilibrium  lengths  of  the  four  compliant  ties.  One  equation  is 
associated  with  the  geometry  constraint  that  must  be  satisfied  in  order  for  the  mechanism 
to  be  assembled.  The  other  three  equations  correspond  to  the  derivative  of  the  potential 
energy  equation  taken  with  respect  to  three  of  the  compliant  tie  lengths. 

Solution  of  these  four  equations  in  order  to  find  sets  of  values  for  the  parameters 
Li  through  L4  that  simultaneously  satisfy  all  four  equations  proved  to  be  a  very  difficult 
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process.  For  this  case,  the  continuity  method  was  used  in  an  attempt  to  find  all 
equilibrium  solutions  for  two  numerical  cases.  The  results  of  the  continuity  method  are 
then  checked  to  first  verify  that  the  set  of  solutions  simultaneously  satisfy  all  four 
equations  and  then  to  detennine  whether  each  solution  set  represents  a  minimum 
potential  energy  and  an  equilibrium  state. 

4.2  Problem  Statement 

The  problem  statement  can  be  explicitly  written  as: 

Given 

a  12,  a34  lengths  of  struts, 

ki,  Loi  spring  constant  and  its  free  length  between  points  4  and  2, 

ki,  L02  spring  constant  and  its  free  length  between  points  3  and  1 , 

k 3 ,  L03  spring  constant  and  its  free  length  between  points  3  and  2, 

k4,  L04  spring  constant  and  its  free  length  between  points  4  and  1 . 

Find  at  equilibrium  position 

Li  length  of  spring  1 , 

L2  length  of  spring  2, 

L3  length  of  spring  3, 

L4  length  of  spring  4, 

It  should  be  noted  that  the  problem  statement  could  be  formulated  in  a  variety  of 
ways.  The  solution  presented  here  uses  L1;  L2,  and  L3  as  the  three  generalized 
parameters,  knowledge  of  which  completely  describes  the  system  (i.e.  the  value  for  L4, 
the  last  remaining  parameter,  can  be  deduced).  The  geometry  equation  that  relates  L4  in 
terms  of  the  three  generalized  parameters  is  derived  first.  Following  this,  a  derivative  of 
the  potential  energy  equation  is  taken  with  respect  to  each  of  the  three  generalized 
parameters  to  yield  an  additional  three  equations. 
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4.3  Development  of  Geometric  and  Potential  Energy  Constraint  Equations 

The  geometry  equation  is  derived  in  a  manner  similar  to  that  presented  in  Section 
3.2.1.  The  parameter  a4i  in  equations  (3.15)  can  be  replaced  by  L4  to  yield  an  equation  in 
the  unknown  terms  Li,  L2,  L3,  and  L4.  This  equation  can  be  rearranged  as 

G!  (Li4  L24)  +  G2  (IV  L22)  +  G3  (L32  U2  L22)  +  G4  (L!2  L22)  +  G5  (L32  L22) 

+  G6  (L22)  +  G7(  L32  Lf  )  +  G8(  U2 )  +  G9(  L34 )  +  Gi0(  L32 )  +  G„  =  0  (4-1) 

where 


G!  =  l 
G2=  1 
G3  =  -l 


^  _  T  2  2  2 

u4  — L4  -  ai2  -  a34 


G5  -  ai22  -  L42 


^  _  T  2  2  2  2 

tj6  -  i-4  a34  -  a34  ai2 


G7  =  -L42  Li2 


n  _  2  2  2  T  2 

u8  -  -  a34  ai2  -  ai2  _l4 

g9=  l42 

_ T  4  2  T  2  2  2  2  T  2 

uio  —  l4  -  ai2  i^4  -  a34  ai2  -  a34  l4 

^_t2  2  2,  4  2  4  2 

uh--l4  a34  ai?  +  a34  ai2  -  aj2  a34  . 


(4-2) 


The  potential  energy  of  the  system  can  be  evaluated  as 

U  =  tk](L]-L„l)2  +  ik,(L,-L0,)2  +  ik3(L3-LM)2+ik4(L4-LM)1  (4-3) 

At  equilibrium,  the  potential  energy  will  be  a  minimum.  This  condition  can  be 
determined  as  the  configuration  of  the  mechanism  whereby  the  derivative  of  the  potential 
energy  taken  with  respect  to  the  lengths  Li,  L2,  and  L3  all  equal  zero,  i.e. 
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^=k,(L,-Ltll)+k4(L4-L(14)^i  =  0  ,  (4-4) 

aLj  aLj 

^J-=k2(L2-Lo!)+k4(L4-Lo4)^±  =  0  ,  (4-5) 

dL2  dL2 

^=kJ(LJ-L02)+k4(L4-L04)^±  =  0  .  (4-6) 

dL3  dL3 


The  derivatives  dL4/dLi,  dL4/dL2,  and  dL4/dL3  are  detennined  by  implicit 
differentiation  of  Equation  4-1  as 


dL4  _ 

T  TT  4  T  2  2  ,  OT  2T  2  T  2T  2  T  2  2  T  2T  2  ,  2T  2  T  2T  2  2  2  ,  T  2  2 

L1[L2  L2  a 34  +2L2  Lj  L2  L3  L2  a12  L2  L4  +a12  L4  L4  L3  a34  a12  +L3  a34  )] 

(4-7) 

dL, 

dL4  _ 

T  (T  2  2  T  2t  2  T  2J  2  2j  2  J  2j  2  2  2  ,  9T  2t  2  ,  t  4  t  2  2  t  2  2 

.L4  ^1^2  ^34  •L'2  ^“'1  ^~“2  ^3  '  ^12  ^1  ^3  ^1  ^34  ^12  ^'*-‘4  ^3  '  ^3  ^3  ^34  ^3  ^12 

T  r  /-»T  2 T  2  2T  2  T  4  ,  x  2X  2  ,  2X  2  x  2  2  ,  T  2X  2  ,  2  2  ,  T  2X  2  x  2  2^-, 

L2  y  2 L2  Lj  +  a34  Lj  Lj  +  L3  Lj  +  aj2  Lj  L4  a34  +  L4  Lj  +  a34  aj2  +  L4  L3  L3  aj2  )J 

(4-8) 

dL2 

dL4  _ 

T  ,T  2  2  y  2  j  2  ,  2,  2  2t  2  y  2y  2  2  2,OT2T2,T4  T2  2  T2  2\ 

■L'4  (L2  ^34  -L'2  •^-'1  J-,2  ^3  ■  ^12  ^1  ^3  •^■'1  ^34  ^12  '  ^-L^4  -L<3  “r  L3  E3  a34  -L/3  ^12  / 

x  r  x  2X  2  x  2X  2  x  2  2  T  2X  2  ,  T  4  iOX  2X  2  x  2  2  2X  2  2  2  2X  2N1 

L3  [  L2  L,  L2  L4  L2  a12  L4  Lj  +L4  +2L4  L3  L4  a34  a12  L4  a34  a12  a34  Lj  )] 

(4-9) 

dL3 

I  /T  2  2  j  2  j  2  T  2T  2  2  j  2  y  2  j  2  2  2  i'll  2T  2  ±1  4  T2  2  T  2  2\ 

*-‘4  V^2  ^34  ^2  ■*“'1  ^2  -^“'3  '  ^12  ^1  ^3  ^1  ^34  ^12  '  ^•*-‘4  -L'3  ^3  ^3  ^34  ^-'3  &12  ) 

Substituting  Equation  4-7  into  Equation  4-4  and  rearranging  gives 
DiCL)  L24)  +  D2(L!3  l22)  +  D3(Lf  L22)  +  D4(L32  l,  l22)  +  D5(L!  l22) 

+  D6(L32  L22)  +  D7(L22)  +  D8(L34)  +  D9(L32)  +  D10  +  Dh(L32  L!3)  +  D^L)3) 

+  D13(L32  L,2)  +  D14(L!2)  +  D15(L34  U)  +  D16(L32  LO  +  D17  (LO  =  0  (4-10) 

where 

D  i  -k4L4  “t-  k  :  1 .  ok  . 

D2  =  -2k4L4  +  2k4Lo4  —  k)L4  , 

D3  =  kiL4Loi  , 

D4  =  -kiL4-k4Lo4  +  k4L4  , 

D5  =  -k4Lo4L4  +  k4L4  +  kjL4a34“  +  k4L4a  1 2  -  k4Lo4ai2  +  k4a34  L4  -  k4Lo4a34" , 

E>6  =  kiL4Loi  , 

D7  =  -kiL4Loia34‘  , 
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D8  -  kiL/tLoi  , 

2  3  2 

D9  =  kiL4Loiai2“-2kiL4  Loi+kiL4Loia34"  , 

2  2 

D  10  -  kiL4Loia34"ai2~  , 

Du  =  kiL4 , 

D12  =  kiL4ai2-  , 

D  13  -  kil^Loi  , 

D14  =  -kiL4Loiai2  , 


D15  -  kiL4 , 

Di6  =  2k  1 L43  +  k4L4^  +  k4Lo4a34-  -  k  1  L4ai22  -  k4a342L4  +  kiL4a342  -  k/iLo/iLzi” 


23  22  22  22  2 

D17  =  -k4ai2  L4  -  kiL4a34~ai2  +  k4L4a34^ai2  +  k4Lo4ai2  L4  -  k4Lo4a34~ai2.  (4-11) 


Substituting  Equation  4-8  into  Equation  4-5  and  rearranging  gives 

E!(Li2L23)  +  E2(L32L23)  +  E3(L23)  +  E4(Li2L22)  +  E5(L32L22)  +  E6(L22)  +  E7(Li4L2) 

+  E8(L32L2)  +  E9(Li2L2)  +  Eio(L34L2)  +  Eh(L32Li2L2)  +  En(U)  +  Ei3(L32) 

+  Ei4(Li2)  +  E15(L34)  +  Ej6  =  0  (4-12) 

where 

E 1  =  -k2L4  —  2k4L4  +  2k4Lo4  , 

E2  =  -k2L4 , 

E3  =  k2L4a342  , 

E4  =  k2L4L02  , 

E5  =  k2L4Lo2  , 

E6  =  -k2L4Lo2a34  , 

E7  =  k4L04  -  k4L4  , 

3  2  3  2  2  2  2 

Eg  =  k4L4  —  k4L04L4-  +  2k2L4  —  !<2L4a34*"  —  k2L4ai2  —  k4L4ai2~  +  k4Lo4ai2"" , 
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E9  —  k4L43  +  kol^a^2  +  k4L4&34~  +  k4L4ai2*"  —  k4Lo4a34_  —  k4Lo4ai22  -  k4Lo4L4~ 

E10  =  k2L4 , 

E 1 1  =  -k2L4  +  k4L4  —  k4Lo4  , 

E12  =  -k4L43a342  +  k4L4a34“ai2~  -  k4Lo4a34"ai22  -  k2L4a34"ai22  +  k4Lo4L42a34"  , 

2  3  2  2 

E13  =  k2L4Lo2ai2  -  2k2L4  L02  +  k2L4Lo2a34~  +  k2L4Lo2Lf 


E14  -  -k2L4Lo2ai2“  , 


E15  =  -k2L4Lo2  , 

Ei6  =  k2L4Lo2a342ai22  .  (4-13) 

Lastly,  substituting  Equation  4-9  into  Equation  4-6  and  rearranging  gives 
FjCLjV)  +  F2(L3L12L22)  +  F3(L33L22)  +  F4(L32L22)  +  F5(L3L22)  +  F6(L22) 

+  F7(L33L12)  +  F8(L3  Lj2)  +  F9(L!2)  +  F10(L35)  +  Fn(L34)  +  F12(L33)  +  F13(L32) 

+  F14(L3)  +  F15  =  0  (4-14) 

where 

F 1  =  k3L4Lo3 , 

F2  =  k4L4  -  k4L04  -  k3L4  , 

F3  =  -k3L4  , 

F4  =  k3L4Lo3  , 

2  2  3  2  2 

F5  =  lc3L4a34'  -  k4L4ai2  +  IC4L4  -  k4Lo4L4  +  k4Lo4ai2‘ , 

F  6  =  -k3L4Lo3a34  , 

F7  =  -k3L4  +  k3L4Lo3  , 

3  2  2  2  2 

Fg  =  k4L4  +  k3L4ai2  -  k4L4a34  -  k4Lo4L4  +  k4Lo4a34  , 

F9  =  -k3L4Lo3ai2  , 

F10  =  k3L4  , 

F 1 1  =  -k3k4Lo3 , 

F12  =  2k3L43  -  l<3L4a34‘  -  k3L4ai2“  +  2k4Lo4L4“  -  2k4L43 , 

3  2  2 

F13  =  -2k3L4  Lq3  +  k3L4Lo3a34”  +  k3L4Lo3ai2-  , 


77 


32  23  4  22  22 

F14  =  k^4  a34~  +  k4an  L4~  +  L4L04L4  +  k3L4a34  a^  +  k4L4a34  ai2- 

2  2  2  2  2  2  5 

-  k4Lo4L4“a34‘'  -  k4Lo4ai2“L4~  -  k4Lo4a34“ai2  -  k4L4~  , 


2  2 

F15  =  k3L4Lo3a34  an"  . 


(4-15) 


4.4  Solution  of  Geometry  and  Energy  Equations 

The  problem  at  hand  is  to  determine  sets  of  values  for  the  parameters  Li,  L2,  L3, 
and  L4  that  will  simultaneously  satisfy  Equations  (4-1),  (4-10),  (4-12),  and  (4-14). 
Sylvester’s  method  is  one  approach  for  solving  this  set  of  equations  for  all  possible 
solutions.  The  parameter  L4  has  been  embedded  in  the  coefficients  of  Equations  (4-1), 
(4-10),  (4-12),  and  (4-14)  and  Table  4-1  shows  the  37  products  of  Li,  L2,  and  L3  that  exist 
in  these  four  equations.  The  Table  indicates  which  products  that  are  multiplied  by  which 
coefficients  for  the  four  equations. 

Table  4-1.  Coefficient  G,  D,  E,  and  F. 
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Table  4-1.  Continued. 


It  is  apparent  from  the  complexity  of  the  four  equations  that  attempting  to  multiply 
the  four  equations  by  different  products  of  Lj,  L2,  and  L3  in  order  to  obtain  a  set  of 
“homogeneous”  equations  where  the  number  of  equations  equals  the  number  of 
unknowns  is  impractical.  For  this  reason,  a  numerical  case  will  be  considered  and  this 
case  will  be  solved  using  the  continuation  approach. 

4.5  Numerical  Case  Study 

The  following  parameters  were  selected  as  a  numerical  example: 

Strut  lengths:  ai2  =  14  in.,  a34  =  12  in. 

Spring  free  lengths  and  spring  constants: 

L01  =  8  in.  ki  =  1  lbf/in. 

L02  =  2.68659245  in.  k2  =  2.0  lbf/in. 

L03  =  3.46513678  in.  k3  =  2.5  lbf/in. 

L04  =  7.3082878  in.  Lj  =  1 .5  lbf/in. 

Equations  (4-1),  (4-10),  (4-12),  and  (4-14)  can  now  be  expressed  respectively  as 

L24Li2  +  (Li4  -  L32Li2  -  440.0  Lj2  +  96.0  L32  -  13824.0)  L22  -  8424.0  Lj2  -  8624.0  L!2  + 
44.0  L32  Li2  +  6773760.0  -  52224.0  L32  +  L34  =  0 


(4-16) 
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-  4.03756830  L24  L,  +  (-1 1520.0  -  5.96243170  L32  Li  +  80.0  Lj2  -  18.07513660.0  Li3  + 
3216.53005200  L,  +  80.0  L32)  L22  +  1 1200.0  L32  +  1960.0  Lj3  -10  L32  Li3  - 
1577.65300520  L32  Li  +  10.0  Li  L34  -  80.0  L34  +  0.2257920  107  -  15680.0  L,2  - 
247420.01098080  L,  +  80.  L32  U2  =  0 

(4-17) 

(-28.07513660  Lf  +  2880.0  -  20.0  L32)  L23  +  (-7737.38625600  +  53.73184900  L32  + 
53.73184900  L!2)  L22  +  (-3187.60655680  L32  +  5696.53005200  Lf  -  15.96243170  L32 
Li2  +  20.0  L34  -  4.03756830  Li4  -  508664.65582080)  L2  +  7522.45886000  L32  + 
0.151652770617600  107  -  53.73184900  L34  -  10531.44240400  Li2  +  53.73184900  L32 

Lj2  =  0 

(4-18) 

-494742.03310080  L3  +  12127.978730000  L32  -  16979.170222000  Lf  +  86.628419500 
L32  Li2  -  86.628419500  L34  -  4307.5 1366000  L33  +  4722.34699480  L3  Lf  +  25.0  L35  - 
25.0  L33  L2  +  (-20.96243170  L3  Lf  +  86.628419500  Lf  +  3212.39344320  L3  -  25.0  L33 
_  12474.492408000  +  86.628419500  L32)  L22  +  0.2445000511968000  107  =  0 

(4-19) 

The  Polynomial  Continuation  Method  was  used  to  find  simultaneous  roots  for  the 
four  equations.  Eighteen  real  solutions  and  490  complex  solutions  were  found.  All  the 
real  and  complex  solutions  were  shown  to  satisfy  Equations  (4-16)  through  (4-19)  with 
neglible  residuals.  The  real  solutions  are  listed  in  Table  4-2. 

Table  4-2.  Eighteen  real  solutions.  (Solutions  that  are  geometrically  impossible  are 


shaded.) 


Lh  in. 

L2,  in. 

L3,  in. 

L4,  in. 

1 

19.028528 

3.624423 

3.108440 

4.306682 

2 

13.000001 

7.999996 

7.016593 

9.999998 

3 

11.632986 

-4.613617 

9.950038 

12.269992 
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Table  4-2.  Continued. 


Li,  in. 

L2,  in. 

L3,  in. 

L4,  in. 

4 

3.141567 

2.156404 

5.241983 

18192258 

5 

-7.985150 

10.049547 

-6.186358 

14.286022 

6 

-0.398136 

4.497517 

10.533965 

15.669134 

7 

-6.868035 

1.105366 

3.486019 

20.971818 

8 

4.200027 

-3.538652 

10.462012 

15.411513 

9 

19.507910 

4.034579 

3.169571 

5.347832 

10 

-16.748403 

8.590494 

2.764444 

2.329743 

11 

4.520731 

0.603385 

-7.518882 

18.852623 

12 

4.077402 

4.708715 

9.316941 

16.646605 

13 

13.473294 

11.823371 

-2.235364 

7.821817 

14 

22.800126 

6.225147 

13.541713 

-5.861200 

15 

14.616270 

-9.580282 

.4.547034 

8.863331 

16 

20.778750 

-0.084697 

4.369638 

6.733159 

17 

-12.774217 

12.496006 

-1.650677 

8.309365 

18 

3.806954 

2.744884 

8.098941 

17.877213 

Cases  1,  4,  7,  9,  10,  11,  14,  16,  and  18  are  not  physically  realizable.  For  example, 


Li  in  case  1  is  greater  than  sum  of  an  and  L4,  and  hence  point  2  cannot  be  constructed. 


The  solutions  in  Table  4-2  can  correspond  to  either  maximum  or  minimum 


potential  energy  configurations.  The  second  derivative  of  the  potential  energy  with 


respect  to  the  variables  Li,  L2,  and  L3  are  now  evaluated  to  identity  the  minimal  potential 


energy  cases. 

Derivatives  of  Equations  (4-4)  through  (4-6)  with  respect  to  Li,  L2,  and  L3, 
respectively  can  be  written  as 


d2U 


dL4  2 


-kj  +k4(—— )“  +k4(L4  -L04) 


d  L, 


2  _JV1  1  4V,  )  1  V  4  04  /  jt2 

1  “Ei 


d  U 
dL, 


,dL 


d2U 


dL 


-k2  +k4(— — )“  +k4(L4  -L04) 
dL, 


2  =k3+k4(^t)!+k4(L4-LM) 


dL 


dL  1 

d2L4 

dL22  ’ 

d2L4 

4V^4  ^04/  1T  t 

dL  3 


(4-20) 
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Table  4-3  shows  the  value  of  the  potential  energy  and  the  second  derivative  for 
each  of  the  remaining  nine  cases.  It  is  seen  that  cases  2,  3,  13,  15,  and  17  have  all 
positive  second  derivative  values  and  represent  possible  equilibrium  solutions.  Each  of 
these  solutions  was  analyzed  to  determine  if  it  truly  corresponds  to  an  equilibrium 
configuration  and  only  cases  2,  3,  13,  and  15  were  indeed  in  equilibrium. 


Table  4-3.  Identification  of  minimum  potential  energy  state.  (Solutions  that  do  not  have 
all  positive  second  derivatives  are  shaded.) 


Case 

Li,  in 

L2,  in 

L3,  in 

L4,  in 

U,  in-lbf 

d2U 
dL 2 

lbf/in 

d2U 

SL2 

lbf/in 

d2U 

dL2 

lbf/in 

2 

13.000 

7.999 

7.016 

9.999 

61.932 

1.146 

3.052 

4.569 

3 

11.632 

-4.613 

9.950 

12.269 

130.923 

0.740 

0.566 

4.065 

5 

-7.985 

10.049 

-6.186 

14.286 

334.931 

-0.470 

1.241 

1.594 

6 

-0.398 

4.497 

10.533 

15.669 

153.432 

3.597 

1.889 

-0.750 

8 

4.200 

-3.538 

10.462 

15.411 

156.415 

0.659 

-16.134 

-2494.334 

12 

4.077 

4.708 

9.316 

16.646 

119.990 

-0.221 

-15.521 

-379.009 

13 

13.473 

11.823 

-2.235 

7.821 

139.276 

4.967 

51.239 

56.665 

15 

14.616 

-9.580 

-4.547 

8.863 

25.420 

2.405 

25.608 

32.573 

17 

-12.774 

12.496 

-1.650 

8.309 

345.474 

10.529 

50.466 

39.862 

Figure  4-2  shows  the  four  feasible  solutions  that  are  in  stable  equilibrium  drawn 
to  scale.  Dashed  lines  are  used  to  identify  the  springs  that  have  negative  spring  lengths. 

It  must  be  noted  that  when  a  spring  has  a  negative  spring  length,  the  two  points  at  the  end 
of  the  spring  want  to  move  towards  one  another  as  the  spring  attempts  to  return  to  its  free 
length. 


Figure  4-2.  Four  equilibrium  solutions 
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A  force  analysis  was  conducted  for  each  of  the  four  cases  and  these  cases  were 
found  to  be  in  static  equilibrium.  Table  4-4  shows  the  calculated  force  in  each  of  the  four 
springs  and  the  two  struts  for  the  equilibrium  cases.  Negative  values  in  the  table 
correspond  to  the  springs  which  have  a  negative  spring  length. 


Table  4-4.  Calculated  forces  in  struts,  and  elastic  ties. 


Case 

Force  in 
Spring 

Lh  lbf 

Force  in 
Spring  L2, 
lbf 

Force  in 
Spring  L3, 
lbf 

Force  in 
Spring  L4, 
lbf 

Force  in 
Strut  ai2, 
lbf 

Force  in 
Strut  a34, 
lbf 

2 

5.000 

10.625 

8.877 

4.036 

11.828 

7.220 

3 

3.632 

-  14.599 

16.212 

7.441 

17.942 

9.253 

13 

5.473 

18.273 

-  14.250 

0.769 

18.5294 

5.8228 

15 

6.616 

-24.533 

-20.030 

2.332 

25.005 

7.6394 

4.6  Conclusions 

The  problem  of  determining  the  equilibrium  configurations  of  a  four-spring  planar 
tensegrity  system  was  formulated  by  developing  four  equations  in  terms  of  the  four 
spring  lengths.  It  was  not  practical  to  use  Sylvester’s  method  to  solve  this  set  of 
equations.  Rather,  the  continuation  method  was  used  for  a  particular  numerical  example 
and  four  equilibrium  configurations  were  obtained.  Cases  may  exist  that  would  yield 
more  than  four  realizable  configurations,  and  this  could  be  investigated  in  future  work. 
The  objective  here,  however,  was  to  determine  if  it  was  tractable  at  all  to  solve  the 
equilibrium  problem  for  multiple  (if  not  all)  solutions,  either  symbolically  or  numerically, 
and  this  goal  was  met. 


CHAPTER  5 

FUTURE  WORK,  SUMMARY,  AND  CONCLUSION 


5.1  Future  Work 
5.1.1  The  3D  Tensegrity  Platform 

As  part  of  future  suggested  research  work,  an  introduction  to  the  position  analysis 
of  a  general  three  dimensional  parallel  platform  device  consisting  of  3  struts,  3  springs, 
and  6  ties  (Figure  5-1)  is  presented  here.  The  problem  statement  is  given  as  follows: 

Given:  length  of  non-compliant  ties,  struts,  and  free  lengths  and  spring  constants 
of  the  compliant  ties. 

Find:  final  length  of  compliant  ties  at  all  equilibrium  configurations. 
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Three  approaches  to  this  problem  are  presented  here: 

1 .  Minimum  potential  energy  analysis. 

2.  Force  balance  analysis. 

3.  Linear  dependence  of  6  connector  lines. 

The  unexpected  complexity,  i.e.  high  degree  polynomials,  encountered  in  the  planar 
analysis  may  indicate  that  the  analysis  of  the  general  three  strut  tensegrity  system  will  be 
very  complex. 

This  is  a  three  degree  of  freedom  system.  One  approach  is  to  define  the  system  by 
six  angles,  subject  to  three  constraints.  That  is,  the  coordinates  of  the  tip  of  each  strut 
(points  4,  5,  and  6)  can  be  defined  by  angles  a;  and  Pi,  i  =  1..3,  which  define  the 
orientation  of  each  of  the  three  struts  relative  to  the  base  triangle.  The  three  constraints 
are  that  the  distances  between  points  4,  5,  and  6  are  given.  The  potential  energy  can  be 
written  in  terms  of  these  six  angles  and  three  equations  can  be  generated  by  equating  the 
partial  derivative  of  the  potential  energy  with  respect  to  three  of  the  angles  to  zero.  These 
three  equations,  together  with  the  three  distance  constraint  equations  would  give  a  set  of 
six  equations  in  the  six  angles. 

In  the  second  approach,  force  balance  equations  at  points  4,  5,  and  6  could  be 
written  as  functions  of  the  lengths  of  the  three  springs.  This  is  complicated  in  that  the 
positions  of  points  4,  5,  and  6  would  have  to  be  determined  from  a  forward  analysis  of 
the  device,  i.e.  the  device  is  equivalent  to  a  3-3  parallel  mechanism,  whose  forward 
analysis  was  solved  by  Griffis  and  Duffy  [27].  The  solution  of  the  forward  analysis  of 
the  3-3  parallel  mechanism  requires  rooting  a  fourth  degree  polynomial  in  the  square  of 
one  parameter,  resulting  in  eight  possible  configurations  for  the  device.  The  nature  of  the 
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forward  analysis  makes  it  infeasible  to  obtain  a  symbolic  solution  to  the  equilibrium 
configuration  problem.  A  numerical  solution  may  be  possible. 

The  third  approach  utilizes  the  fact  that  if  the  device  is  in  equilibrium  with  no 
external  forces  applied,  then  the  Pliicker  coordinates  of  the  lines  along  the  three  struts  and 
the  three  springs  must  be  linearly  dependent.  For  this  solution,  it  is  then  necessary  to 
obtain  the  Pliicker  coordinates  of  these  lines  as  functions  of  the  three  spring  lengths.  This 
again  requires  a  forward  analysis  of  the  3-3  parallel  mechanism.  Further,  only  a  single 
equation  results  from  equating  the  determinant  of  the  6X6  matrix  of  line  coordinates  to 
zero.  This  leads  one  to  believe  that  the  fact  that  the  lines  must  be  linearly  dependent  is  a 
necessary,  but  not  sufficient  condition  for  equilibrium.  Further  analysis  is  necessary. 

5,1.2  The  2D  Tensegrity  structures 

Cases  of  two  dimensional  tensegrity  structures  proved  to  be  mathematically 
challenging  in  spite  of  their  simple  configurations.  Different  methods  were  tried  to  solve 
these  problems.  It  was  expected  the  final  closed  form  solution  be  a  simple  mathematical 
expression.  In  each  case  the  intermediate  mathematical  expressions  become  complicated 
to  solve  and  the  relevant  matrices  become  large  to  manipulate.  Other  mathematical 
procedures  might  exist  to  solve  these  problems  with  less  complexity. 

5.2:  Summary  and  Conclusion 

In  this  research,  all  possible  combinations  of  two  dimensional  tensegrity  structures 
consisting  of  two  struts,  springs,  and  non-compliant  members  were  analyzed.  A  closed 
form  solution  could  only  be  obtained  for  case  of  2  spring,  2  struts,  and  2  non-compliant 
members.  Finding  a  simple  solution  approach  for  a  closed  form  solution  to  the  other  two 
dimensional  tensegrity  structures  consisting  of  either  3  or  4  springs  was  not  possible  due 
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to  the  complexity  of  mathematical  solutions.  A  numerical  solution  using  the  Continuity 
Method  was  used  to  solve  individual  numerical  examples. 

The  contribution  of  this  research  was  to  set  a  basic  mathematical  foundation  for  the 
closed  fonn  solutions  of  tensegrity  structures  and  also  to  present  numerical  approaches 
where  a  closed  form  solution  could  not  be  obtained.  The  contribution  of  this  research 
leads  to  a  better  understanding  of  the  complexity  of  these  devices  and  hopes  to  lay  a 
ground  work  for  future  analyses  of  more  complex  systems. 


APPENDIX 
SAMPLE  OUTPUT 


The  following  coefficients  belong  to  the  numerical  example  of  a  Maple  sessions 
programming. 

(Pl  Cl2  +  p2  Cl  +  p3)  C43  +  (P4  Cl3  +  p5  Cl2  +  p6  Cl  +  p7)  C42  +  (P8  Cl4  +  p9  Cl3  +  pio  Cl2  + 

Pi  1  Cl  +  P12)  C4  +  (p  1 3  Cl4  +  P14  Cl3  +  p  1 5  Cl2  +  Pi6  Cl  +  P17)  =  0  (3-65) 

(qi  ci  +  q2)  c44  +  (q3  Ci2  +  q4  Ci  +  q5)  c43  +  (q6  Ci3+  q7  Ci2  +  q8  Cj  +  q9)  c42  +  (qio  cf3  +  qn 
ci2  +  qi2  ci  +  qi3)  c4+  (qi4  cf’  +  qi5  Ci2  +  qi6  Ci  +  qn)  =  0  (3-66) 

(ri  ci  +  r2)  c42  +  (r3  Ci2  +  r4  Ci  +  r5)  c4  +  (r6  Ci2  +  r7  ci  +  r8)  =  0  (3-69) 


># - 

>  #  Numerical  example 

al2  :=  14  a34  :=  12  a41  :=  10 
L01  :=  8  kl  :=  1 
L02:=  2  k2  :=  2.687 
L03:=  2.5  k3  :=  3.465 


pl  :=  0.5773663961  10 15  L32  -  0.2080248063  1016 13  +  0.1873780685  1016 

p2  :=  0.1110100300  1016  L32  -  0.4198962128  1016  L3  +  0.3961707733  1016 

p3  :=0.5335964670  1015  L32  -  0.21 14119475  10 16  L3  +  0.2094045517  1016 

p4  :=  0.4853912146  10 15  L32  -  0.2087909024  1016 13  +  0.2186077465  1016 
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p5  :=  -0.1733540051  1013  L34  +  0.7456817937  10 13  L33  +  0.2423387462  1016Z32 
-  0.9462608682  1016 13  +  0.9212755033  1016 

p6  :=  -0.3864997467  1013  L34  +  0.1452389258  10 14  L33  +  0.3406073884  1016Z32 
-0.1266504024  1017  L3  +  0.1 176377359  1017 

p7  :=  -0.21 13487353  1013  L34  +  0.7162795261  1013  ZJ3  +  0.1466145513  1016ZJ2 
-0.5300984714  1016 13  +  0.4722921577  1016 

p8  :=  0.1020169484  10 15  L32  -  0.5100847420  10 15  L3  +  0.6376059274  1015 

p9  := -0.7286924884  1012  L34  +  0.3643462442  1013  ZJ3  +  0.1068305616  1016Z32 
-0.4737062676  1016 13  +  0.5137282043  1016 

plO  :=  0.1301236586  1010  ZJ6  -  0.6506182936  10 10  L35  -  0.4747632005  1013Z34 
+  0.2052962276  1014  L33  +  0.3082075504  1016  L32  -  0.1236996613  1017L3 
+  0.1219206631  1017 

pll  :=  0.3300435984  1010  L36  -  0.1289836349  1011  L35  -  0.7989672476  1013Z34 
+  0.2941695331  1014  L33  +  0.3455740168  1016  L32  -  0.1248534533  1017Z3 
+  0.1113665191  1017 

pl2  :=  0.2092793462  1010  L36  -  0.5893635686  1010  L35  -  0.3983027155  1013U4 
+  0.1246857458  1014  L33  +  0.1339547765  1016  L32  -  0.4345164946  1016ZJ 
+  0.3440522646  1016 

p!3  :=  0.1037172309  1015  L32  -  0.5185861544  10 15  L3  +  0.6482326928  1015 

p!4  := -0.7408373632  1012  L34  +  0.3704186816  1013  L33  +  0.5844049572  1015ZJ2 
-0.2657927757  1016  L3  +  0.2963349453  1016 

pl5  :=  0.1322923863  1010  L36  -  0.6614619316  1010Z35  -  0.3034952946  1013Z34 
+  0.1316433327  1014  Z33  +  0.1235296500  1016  L32  -  0.4981478398  1016L3 
+  0.4841901337  1016 

pl6  :=  0.3355443251  1010  L36  -  0.131 1333621  1011  L35  -  0.4128749122  1013Z34 
+  0.1489517912  1014  L33  +  0.1 159476235  1016  Z32  -  0.4015136920  1016ZJ 
+  0.3326208569  1016 

^77:=  0.2121019430  1010  L36  -  0.5991862948  1010  ZJ5  -  0.1862278592  1013Zi4 
+  0.5272839394  1013  L33  +  0.4069170754  1015  Zi2  -  0.1 160024667  1016ZJ 
+  0.8167029969  1015 
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ql  :=  0.8744309862  1014  L32  -  0.4372154932  1015 13  +  0.5465193664  1015 
q2  \=  0.9243984712  1014  Z32  -  0.4621992356  10 15 13  +  0.5777490445  1015 
g3  1=0.4669520350  1015  Z32  -  0.2187549572  1016Z3  +  0.2550423709  1016 

q4  := -0.7286924880  1 0 12  Z34  +  0.3643462442  10 13  L33  +  0.9975040455  1015Z32 
-0.4749518798  1016  L3  +  0.5610932160  1016 

q5  := -0.7703320590  1012  Z34  +  0.3851660295  1013  L33  +  0.5326615535  1015Z32 
-0.2576221616  1016  L3  +  0.3081328237  1016 

q6  :=  0.6233888280  1015  L32  -  0.2723886871  10 16  L3  +  0.2975494328  1016 

q7  :=  -0.1945633479  1013  Z34  +  0.91 14789880  10 13  Z33  +  0.2467565040  1016Z32 
-0.1115537420  1017  L3  +  0.1253958324  1017 

q8  :=  0.1518109351  10 10  Z36  -  0.7590546756  1 0 10  L35  -  0.4948249499  1013Z34 
+  0.2326400716  1014  Z33  +  0.2978451536  1016  Z32  -  0.1366694272  1017Z3 
+  0.1548714436  1017 

q9  :=  0.1604858457  10 10  Z36  -  0.8024292286  1010  L35  -  0.3056661932  1013Z34 
+  0.1440713627  1014  ZJ3  +  0.1 127231047  1016  L32  -  0.5199241 185  1016Z3 
+  0.5873781950  1016 

qlO  :=  0.1238733084  1016  L32  -  0.5475596549  1016  ZJ  +  0.6050171798  1016 

qll  := -0.4022366239  1013  L34  +  0.1777220459  1014 L33  +  0.35541 16422  1016Z32 
-0.1568854809  1017 13  +  0.1720685863  1017 

ql2  :=  0.3260393580  1010  L36  -  0.1447643958  1011  L35  -  0.7932329962  1013Z34 
+  0.3478741954  1014Z33  +  0.3352073482  1016 L32  -  0.1470298991  1017Z3 
+  0.1588185279  1017 

q!3  :=  0.3446701784  1010  L36  -  0.1530366470  1011  L35  -  0.3895030685  1013Z34 
+  0.1691392672  1014  L33  +  0.1035832980  1016  L32  -  0.4479338270  1016L3 
+  0.4707584806  1016 

q!4  :=  0.6153702088  1015  Z32  -  0.2751416753  1016 13  +  0.3075503998  1016 

^75  :=  -0.2075805867  1013  Z34  +  0.8667876209  1013  L33  +  0.1553312052  1016Z32 
-0.6722506098  1016Z3  +  0.7212669096  1016 

ql6  :=  0.1750560045  1010  L36  -  0.6792486624  1010  L35  -  0.3717087746  1013Z34 
+  0.1514057165  1014  Z33  +  0.1284193314  1016  L32  -  0.53468261 16  1016Z3 
+  0.5463422536  1016 
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ql7  := -0.7180628716  1010  L35  -  0.1603893518  1013  L34  +  0.6318953271  1013T33 
+  0.1830521647  1010  L36  -  0.1390169719  1016 13  +  0.1348526898  1016 
+  0.3482593233  1015L32 


rl  \=  8  a!2  a41  a34 2 

r2  :=  4  a412  a342  +  4  a342  a!2 2 

r3  :=  8  a34  a47  a722 

:=  12  al2  a41 2  a34  -  4  al2  a34  L3 2  +  4  al2  a34 3  +  4  al23  a34 
r5  :=  -4  a34  a41  L32  +  4  a413  a34  +  4  a34  a41  al2 2  +  4  a343  a41 
r6  :=  4  a4i2  al2 2  +  4  a342  al2 2 

r7  :=  4  ai23  a47  -  4  ai2  a47  Z32  +  4  oi2  a41  a34 2  +  4  ai2  a413 

r8  \=  L3 4  +  (-2  a342  -2a412-2  al22)L3 2  +  a414  +  a344  +  al24  +  2a412  a34 2 
+  2a412  a!22  -2a342  a!22 


The  following  portion  of  Maple  program  presents  the  three  equations  used  in 
continuation  method. 

(3-71)  =  eqnl_for_cont 
(3-72)  =  eqn3_16_sq 
(3-73)  =  eqn3_17_sq 


>  #  3  springs  cosine  of  thetal,  cosine  of  theta2  and  L23,  Continuation  method 


>  a41 : =  10 


al2 : =  14.0 


a34 : =  12  ; 

a41  :=  10 

al2  :=  14.0 
a34  :=  12 


>  k2  4 : =  1  ;  k31:=  2.0  ;  k23:=  2.5  ; 

k24  :=  1 


k31  :=  2.0 
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k23  :=  2.5 


> L024 : =  8.0;  L031:=  2.686592430;  L023:=  3.465136806; 

L024  :=  8.0 
L031  :=  2.686592430 
L023  :=  3.465136806 


>  eqnl_f or_cont : =  eqnl_f or_cont  ; 

eqnlJor_cont  :=  61600.000  cl  +  42624.00  c4 2  -  168.0  cl  c4  L23 2  +  ^L234 

-  220.0000000  L232  +  47824.00  cl2  +  52800.00  c4  -  140.0  cl  L232 

+  47040.00  cl2c4  +  40320.0  cl  c42  +  107520.000  cl  c4  +  20176.00000 

-  120  c4  L23 2 

> 

>  eqn3_l  6_sq :  =  eqn3_l  6_sq  ; 

eqn3_16_sq  :=  0.2094210873  10 16  c43  +  0.1119744000  10 10  L236  c4 

+  0.1131754495  1010  L236  -  0.4383037685  1010  L235  +  0.3326471226  1016ci 
+  0.4723294525  1016  c42  +  0.5399138304  1014  L232  cl4 
+  0.1817487566  1016ci  c4  L231  -  0.3619570961  10 16  L23  cl2 
+  0.3038915898  1015  L232  cl3  -  0.4221575075  10 13  L234  cl  c4 

-  0.9917251153  1012  L234  +  0.2153950907  1015  L232  +  0.4842283683  1016c72 

-  0.8485560960  1015  L23  +  0.3440794329  1016c4 

+  0.1488688451  1014  L233  cl2  c4  -  0.2488945515  1013  L234  cl2  c4 
+  0.6376562766  1015  ci4  c4  +  0.1275251673  1016  L231  cl2  c42 
+  0.1614293490  1 0^ 16  L232  cl 2  c4  +  0.1795299633  10 16  L232  c42  cl 

-  0.3793305600  1012  L234  cl3  c4  -  0.4694401018  10 10  L235  cl2  c4 
+  0.6773760000  109  L236  cl2  c4  -  0.2043740160  1013  L23 4  cl  c42 

-  0.9388802036  10 10  L235  cl  c4  -  0.1535416328  10 16  L23  c43 

+  0.1741824000  10 10  L236  cl  c4  +  0.5580683989  1015  L231  cl3  c4 
+  0.2141076067  1 0^ 14 L233  cl  c4  +  0.2628864571  1013  L233  cl3  c4 
+  0.9213482524  1016  cl2  c42  +  0.2186250090  1016  cl 3  c42 
+  0.6089439462  10 15  ci  L232  -  0.9103933440  10' 12  L234  cl 2  c42 
+  0.5137687713  1016  cl2  c4  +  0.1055837852  1014  L233  cl  c42 
+  0.6445019392  1015  L232  cl2  +  0.5407949973  1013  L233  cl2  c42 
+  0.3857073162  1013  L233  +  0.1873928649  10 16  cl2  c43 
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+  0.3962020571  1016  cl  c43  -  0.3741750572  10 15  L23  cl4 
+  0.5310627840  1014  L232  cl4  c4  -  0.431 1 184608  10 10  L235  c4 
+  0.2963583455  1016  cl3  +  0.2814309827  10 15  L232  c43 
+  0.2549101363  1015  L231  cl3  c42  +  0.1219302908  1017ci2c4 
+  0.5868135383  1015  L232  c43  cl  -  0.21 19474277  10 13  L234  c4 

-  0.3053174043  1016  L23  c43  cl  -  0.3856527360  10 12  L234  cl3 

-  0.9545282070  1010  L235  cl  -  0.1 122729984  1013  L234  c42 
+  0.6482838809  1015  cl4  +  0.9545282070  1013  L233  cl2 

+  0.7738436783  1015  L232  c42  +  0.1 176470253  10 11  cl  c42 

-  0.4772641035  1010  L235  cl2  +  0.1 1 13753133  10 11  cl  c4 
+  0.6886656000  109Z236  cl2  +  0.1085434933  10 14  L233  cl 

-  0.1924328866  1016L23  cl3  -  0.2180321636  10 13  L234  cl 
+  0.5224006096  1013  L233  c42  -  0.2927947095  1016Z23ci 

-  0.1514225992  1016L23  cl2  c43  +  0.1770854400  10 10  L236  cl 

-  0.6875307066  1016  L23  cl2  c42  -  0.9212468531  10 16  L23  cl  c42 
+  0.3058921636  1015  L232  cl2  c43  -  0.3171517401  10 16  L23  c4 

-  0.3859569282  1016  L23  c42  +  0.9104915321  10 13  L233  c4 

-  0.9090163022  10 16  L23  cl  c4  +  0.8167674885  1015 

-  0.1589435265  1013  L234  cl2  +  0.7079662564  10 15  c4  L232 

-  0.3432245583  1016  L23  cl3  c4  -  0.8985008439  10 16  L23  cl2  c4 
+  0.2672678980  1013  L233  cl 3  -  0.3680410399  1015  L23  cl4  c4 

-  0.1514225992  1016L23  cl3  c42 
> eqn3_17_sq:=  eqn3_17_sq  ; 

eqn3_l  7 _sq  :=  0.3081571556  1016  c43  +  0.1949337600  10 10  L236  c4 

+  0.1117043200  1010L236- 0.5628934554  1010  L235  +  0.5463853959  10 16  cl 
+  0.5874245780  1016  c42  +  0.191 1253794  10 16  cl  c4  L232 

-  0.5115729560  1016L23  cl 2  +  0.3470298604  10 15  L232  cl3 

-  0.4526879307  10 13  L234  cl  c4  -  0.9760319345  1012Z234 

+  0.2101294121  1015  L232  +  0.7213238649  1016  cl2  -  0.1089761729  1016Z23 
+  0.4707956544  1016  c4  +  0.1339260450  10 14  L233  cl2  c4 

-  0.2286520320  1013  L234  cl2  c4  +  0.1380473656  1016  L23 2  cl2  c42 
+  0.2019950461  1016  L232  cl2  c4  +  0.1648149759  10 16 L232  c42  cl 

-  0.2693873525  1013  L234  cl  c42  -  0.1095360238  1011  L235  cl  c4 

-  0.1889778074  1016L23  c43  +  0.1843968000  10 10  L236  cl  c4 

+  0.7038352097  1 0^ 15  L232  cl 3  c4  +  0.2635290684  10 14  L233  cl  c4 
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+  0.1254057344  1017  cl 2  c4 2  +  0.2975729290  1016  cl 3  c42 
+  0.7563577221  1015  cl  L232  -  0.1062125568  1013  Z234  ci2  c42 
+  0.6050649557  1016  cl 3  c4  +  0.1721030006  1014Z233  cl  c4 2 
+  0.8939345721  1 0 15  L232  ci 2  +  0.6747419064  1013  Z233  ci2  c42 
+  0.4953462407  1013  L233  +  0.2550625106  1016  cl2  c43 

+  0.5611375232  1016  cl  c43  -  0.1 157952251  1011  Z235  c4  +  0.3075746857  1016  ci3 

-  0.3334902483  1015  c44  Z23  +  0.5465625225  1015  c44ci 

+  0.2779085403  1 0 13  c43  Z233  +  0.2839097906  1015  L232 c43 
+  0.3568741909  1015  L232  ci3  c42  +  0.1720821738  1017  cl2  c4 
+  0.5380397453  1015  L232  c43  cl  +  0.1348633385  1016 

-  0.2234879229  10 13  L234  c4  -  0.3499544516  10 16  L23  c43  cl 
+  0.5777946668  1015  c44  -  0.5324667821  10 10  L235  cl 

-  0.4010065920  1012  L234  c43  +  0.2628864571  10 13  c43  L233  cl 

-  0.1660829642  1013  L234  c42  -  0.3154637484  10 15  c44  L23  cl 
+  0.6662224379  1013  L233  cl2  +  0.6209287703  1015  L232  c42 

+  0.1548836731  10 11  cl  c42  +  0.1588310692  1017ci  c4 
+  0.8354304000  109  L236  c42  -  0.5789761256  10 10  L235  c42 
+  0.1172863060  1014  L233  cl  -  0.2066277815  10 16  L23  cl3 

-  0.2220502329  1013  L234  ci  +  0.1065316071  1014  L233  c42 

-  0.4129741774  1016Z23  cl  -  0.1619380575  10 16  L23  cl2  c43 
+  0.1075648000  1010  L236  cl  -  0.8361892426  10 16  L23  cl2  c42 

-  0.1021681926  10 11 L23  cl  c42  +  0.2549101363  1015  Z232  ci2  c43 

-  0.3430626536  1016  L23  c4  -  0.3890719565  1016  L23  c42 

+  0.7902720000  109  L236  c42  cl  -  0.5476801 188  1010  L235  c42  cl 
+  0.1289186839  1014  L233  c4  +  0.4812079104  10 14  L231  c44 

-  0.3793305600  1012  L234  c43  cl  +  0.4551966720  10 14  L232  c44  cl 

-  0.1116496308  10 11 L23  cl  c4  -  0.1221936128  10 13  L234  cl2 
+  0.5952551997  1015  c4  L232  -  0.4127784728  10 16  L23  cl3  c4 

-  0.1185533798  1017  L23  cl2  c4  -  0.2061029824  10 16 L23  cl3  c42 

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

>  #  end 
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